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Abstract 

This is a review of neutrino astronomy anchored to the observational fact that Nature accelerates 
protons and photons to energies in excess of 10 20 and 10 13 eV, respectively. 
| Although the discovery of cosmic rays dates back close to a century, we do not know how and 
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where they are accelerated. There is evidence that the highest energy cosmic rays are extra-galactic 
- they cannot be contained by our galaxy's magnetic field anyway because their gyroradius far 



, exceeds its dimension. Elementary elementary-particle physics dictates a universal upper limit 

> ■ 

t**- ' on their energy of 5 x 10 19 eV, the so-called Greisen-Kuzmin-Zatsepin cutoff; however, particles 

in excess of this energy have been observed by all experiments, adding one more puzzle to the 

speculations about the sources. 

There is a realistic hope that the oldest problem in astronomy will be resolved soon by ambitious 
experimentation: air shower arrays of 10 4 km 2 area, arrays of air Cerenkov detectors and, the 



j§ ! subject of this review, kilometer-scale neutrino observatories. 

We will review why cosmic accelerators are also expected to be cosmic beam dumps producing 
associated high-energy photon and neutrino beams. We will work in detail through an example 
of a cosmic beam dump, gamma ray bursts. These are expected to produce neutrinos from MeV 
to EeV energy by a variety of mechanisms. We will also discuss active galaxies and GUT-scale 
remnants, two other classes of sources speculated to be associated with the highest energy cosmic 
rays. Gamma ray bursts and active galaxies are also the sources of the highest energy gamma rays, 
with emission observed up to 20TeV, possibly higher. 

The important conclusion is that, independently of the specific blueprint of the source, it takes 
a kilometer-scale neutrino observatory to detect the neutrino beam associated with the highest 
energy cosmic rays and gamma rays. We also briefly review the ongoing efforts to commission such 
instrumentation. 
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I. THE HIGHEST ENERGY PARTICLES: COSMIC RAYS, PHOTONS AND 
NEUTRINOS 



A. The New Astronomy 

Conventional astronomy spans 60 octaves in photon frequency, from 10 4 cm radio- waves 
to 10~ 14 cm gamma rays of GeV energy; see Fig. 1. This is an amazing expansion of the power 
of our eyes which scan the sky over less than a single octave just above 10~ 5 cm wavelength. 
This new astronomy probes the Universe with new wavelengths, smaller than 10 -14 cm, or 
photon energies larger than 10 GeV. Besides the traditional signals of astronomy, gamma 
rays, gravitational waves, neutrinos and very high-energy protons become astronomical mes- 
sengers from the Universe. As exemplified time and again, the development of novel ways 
of looking into space invariably results in the discovery of unanticipated phenomena. As is 
the case with new accelerators, observing only the predicted will be slightly disappointing. 



FIG. 1: The diffuse flux of photons in the Universe, from radio waves to GeV-photons. Above tens 
of GeV, only limits are reported although individual sources emitting TeV gamma rays have been 
identified. Above GeV energy, cosmic rays dominate the spectrum. 
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Why pursue high-energy astronomy with neutrinos or protons despite considerable instru- 
mental challenges? A mundane reason is that the Universe is not transparent to photons 
of TeV energy and above (units are: GeV/TeV/PeV/EeV/ZeV in ascending factors of 10 3 ). 
For instance, a PeV energy photon cannot deliver information from a source at the edge of 
our own galaxy because it will annihilate into an electron pair in an encounter with a 2.7 
Kelvin microwave photon before reaching our telescope. In general, energetic photons are 
absorbed on background light by pair production 7 + 7bk g nd — > e + + e~ of electrons above 
a threshold E given by 

4£e~(2m e ) 2 , (1) 

where E and e are the energy of the high-energy and background photon, respectively. 
Eq. (1) implies that TeV-photons are absorbed on infrared light, PeV photons on the cosmic 
microwave background and EeV photons on radio- waves; see Fig. 1. Only neutrinos can 
reach us without attenuation from the edge of the Universe. 

At EeV energies, proton astronomy may be possible. Near 50 EeV and above, the arrival 
directions of electrically charged cosmic rays are no longer scrambled by the ambient mag- 
netic field of our own galaxy. They point back to their sources with an accuracy determined 
by their gyroradius in the intergalactic magnetic field B: 

d dB 

6 = = — , (2) 

BE 
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where d is the distance to the source. Scaled to units relevant to the problem, 
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Speculations on the strength of the inter-galactic magnetic field range from 1CT 7 to 
1CT 12 Gauss in the local cluster. For a distance of 100 Mpc, the resolution may therefore 
be anywhere from sub-degree to nonexistent. It is still possible that the arrival directions 
of the highest energy cosmic rays provide information on the location of their sources. Pro- 
ton astronomy should be possible; it may also provide indirect information on intergalactic 
magnetic fields. Determining the strength of intergalactic magnetic fields by conventional 
astronomical means has been challenging. 
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B. The Highest Energy Cosmic Rays: Facts 

In October 1991, the Fly's Eye cosmic ray detector recorded an event of energy 3.0 ±(j; 54 
xl0 20 eV 0. This event, together with an event recorded by the Yakutsk air shower array 
in May 1989 0, of estimated energy ~ 2 x 10 20 eV, constituted (at the time) the two 
highest energy cosmic rays ever seen. Their energy corresponds to a center of mass energy 
of the order of 700 TeV or ~ 50 Joules, almost 50 times the energy of the Large Hadron 
Collider (LHC). In fact, all active experiments |5] have detected cosmic rays in the vicinity of 
100 EeV since their initial discovery by the Haverah Park air shower array M . The AGAS A 
air shower array in Japan 0] has now accumulated an impressive 10 events with energy in 
excess of 10 20 eV §]. 

The accuracy of the energy resolution of these experiments is a critical issue. With a 
particle flux of order 1 event per km 2 per century, these events are studied by using the 
earth's atmosphere as a particle detector. The experimental signature of an extremely high- 
energy cosmic particle is a shower initiated by the particle. The primary particle creates 
an electromagnetic and hadronic cascade. The electromagnetic shower grows to a shower 
maximum, and is subsequently absorbed by the atmosphere. 

The shower can be observed by: i) sampling the electromagnetic and hadronic components 
when they reach the ground with an array of particle detectors such as scintillators, ii) 
detecting the fluorescent light emitted by atmospheric nitrogen excited by the passage of 
the shower particles, iii) detecting the Cerenkov light emitted by the large number of particles 
at shower maximum, and iv) detecting muons and neutrinos underground. 

The bottom line on energy measurement is that, at this time, several experiments using 
the first two techniques agree on the energy of EeV-showers within a typical resolution of 
25%. Additionally, there is a systematic error of order 10% associated with the modeling of 
the showers. All techniques are indeed subject to the ambiguity of particle simulations that 
involve physics beyond the LHC. If the final outcome turns out to be an erroneous inference 
of the energy of the shower because of new physics associated with particle interactions at 
the Aq CD scale, we will be happy to contemplate this discovery instead. 

Could the error in the energy measurement be significantly larger than 25%? The answer 
is almost certainly negative. A variety of techniques have been developed to overcome the 
fact that conventional air shower arrays do calorimetry by sampling at a single depth. They 
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FIG. 2: The cosmic ray spectrum peaks in the vicinity of 1 GeV and has features near 10 15 and 
10 19 eV referred to as the "knee" and "ankle" in the spectrum, respectively. Shown is the flux of 
the highest energy cosmic rays near and beyond the ankle measured by the AGASA experiment. 
Note that the flux is multiplied by E 3 . 



also give results within the range already mentioned. So do the fluorescence experiments 
that embody continuous sampling calorimetry. The latter are subject to understanding the 
transmission of fluorescent light in the dark night atmosphere — a challenging problem 
given its variation with weather. Stereo fluorescence detectors will eventually eliminate 
this last hurdle by doing two redundant measurements of the same shower from different 
locations. The HiRes collaborators have one year of data on tape which should allow them 
to settle energy calibration once and for all. 

The premier experiments, HiRes and AGASA, agree that cosmic rays with energy in 
excess of 10 EeV are not galactic in origin and that their spectrum extends beyond 100 EeV. 
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FIG. 3: As in Fig. 2, but as measured by the HiRes experiment. 

They disagree on almost everything else. The AGASA experiment claims evidence that 
the highest energy cosmic rays come from point sources, and that they are mostly heavy 
nuclei. The HiRes data do not support this. Because of such low statistics, interpreting the 
measured fluxes as a function of energy is like reading tea leaves; one cannot help however 
reading different messages in the spectra (see Fig. 2 and Fig. 3). 



C. The Highest Energy Cosmic Rays: Fancy 

1. Acceleration to > 100 EeV? 

It is sensible to assume that, in order to accelerate a proton to energy E in a magnetic 
field B, the size R of the accelerator must be larger than the gyroradius of the particle: 

R > R g yro = | ■ (4) 

That is, the accelerating magnetic field must contain the particle orbit. This condition yields 
a maximum energy 

E = yBR (5) 
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TABLE I: Requirements to generate the highest energy cosmic rays in astrophysical sources. 
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GRB 7 > 10 2 B = 


10 12 G M = 


M sun 



by dimensional analysis and nothing more. The 7-factor has been included to allow for the 
possibility that we may not be at rest in the frame of the cosmic accelerator. The result 
would be the observation of boosted particle energies. Theorists' imagination regarding 
the accelerators has been limited to dense regions where exceptional gravitational forces 
create relativistic particle flows: the dense cores of exploding stars, inflows on supermassive 
black holes at the centers of active galaxies, annihilating black holes or neutron stars. All 
speculations involve collapsed objects and we can therefore replace R by the Schwartzschild 
radius 

R ~ GM/c 2 (6) 

to obtain 

E oc -fBM. (7) 

Given the microgauss magnetic field of our galaxy, no structures are large or massive enough 
to reach the energies of the highest energy cosmic rays. Dimensional analysis therefore limits 
their sources to extragalactic objects; a few common speculations are listed in Table 1. 

Nearby active galactic nuclei, distant by ~ 100 Mpc and powered by a billion solar mass 
black holes, are candidates. With kilogauss fields, we reach 100 EeV. The jets (blazars) 
emitted by the central black hole could reach similar energies in accelerating substructures 
(blobs) boosted in our direction by Lorentz factors of 10 or possibly higher. The neutron 
star or black hole remnant of a collapsing supermassive star could support magnetic fields 
of 10 12 Gauss, possibly larger. Highly relativistic shocks with 7 > 10 2 emanating from the 
collapsed black hole could be the origin of gamma ray bursts and, possibly, the source of 
the highest energy cosmic rays. 
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The above speculations are reinforced by the fact that the sources listed are also the 
sources of the highest energy gamma rays observed. At this point, however, a reality check 
is in order. The above dimensional analysis applies to the Fermilab accelerator: 10 kilogauss 
fields over several kilometers corresponds to 1 TeV. The argument holds because, with opti- 
mized design and perfect alignment of magnets, the accelerator reaches efficiencies matching 
the dimensional limit. It is highly questionable that nature can achieve this feat. Theorists 
can imagine acceleration in shocks with an efficiency of perhaps 10%. 

The astrophysics problem of obtaining such high-energy particles is so daunting that many 
believe that cosmic rays are not the beams of cosmic accelerators but the decay products 
of remnants from the early Universe, such as topological defects associated with a Grand 
Unified Theory (GUT) phase transition. 

2. Are Cosmic Rays Really Protons: the GZK Cutoff? 

All experimental signatures agree on the particle nature of the cosmic rays — they look 
like protons or, possibly, nuclei. We mentioned at the beginning of this article that the 
Universe is opaque to photons with energy in excess of tens of TeV because they annihilate 
into electron pairs in interactions with the cosmic microwave background. Protons also 
interact with background light, predominantly by photoproduction of the A-resonance, i.e. 
V + Icmb — > A — > 7r + p above a threshold energy E p of about 50 EeV given by: 

2E p e > (ml - m 2 p ) . (8) 

The major source of proton energy loss is photoproduction of pions on a target of cosmic 
microwave photons of energy e. The Universe is, therefore, also opaque to the highest energy 
cosmic rays, with an absorption length of 

A 7P = (n C MB 0p +7CMB ) _1 (9) 
= lOMpc, (10) 

when their energy exceeds 50 EeV. This so-called GZK cutoff establishes a universal upper 
limit on the energy of the cosmic rays. The cutoff is robust, depending only on two known 
numbers: u C mb = 400 cm" 3 and cr p+7CMB = 10~ 28 cm 2 |f| P, |i~0l, pTT| . 

Protons with energy in excess of 100 EeV, emitted in distant quasars and gamma ray 
bursts, will lose their energy to pions before reaching our detectors. They have, nevertheless, 
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been observed, as we have previously discussed. They do not point to any sources within 
the GZK-horizon however, i.e. to sources in our local cluster of galaxies. There are three 
possible resolutions: i) the protons are accelerated in nearby sources, ii) they do reach us 
from distant sources which accelerate them to even higher energies than we observe, thus 
exacerbating the acceleration problem, or iii) the highest energy cosmic rays are not protons. 

The first possibility raises the challenge of finding an appropriate accelerator by confining 
these already unimaginable sources to our local galactic cluster. It is not impossible that all 
cosmic rays are produced by the active galaxy M87, or by a nearby gamma ray burst which 
exploded a few hundred years ago. 



Stecker [12| has speculated that the highest energy cosmic rays are Fe nuclei with a 
delayed GZK cutoff. The details are complicated but the relevant quantity in the problem 
is 7 = E/AM, where A is the atomic number and M the nucleon mass. For a fixed observed 
energy, the smallest boost above GZK threshold is associated with the largest atomic mass, 
i.e. Fe. 



3. Could Cosmic Rays be Photons or Neutrinos? 

Topological defects predict that the highest energy cosmic rays are predominantly pho- 
tons. A topological defect will suffer a chain decay into GUT particles X and Y, that 
subsequently decay to familiar weak bosons, leptons and quark or gluon jets. Cosmic rays 
are, therefore, predominately the fragmentation products of these jets. We know from accel- 
erator studies that, among the fragmentation products of jets, neutral pions (decaying into 
photons) dominate, in number, protons by close to two orders of magnitude. Therefore, if 
the decay of topological defects is the source of the highest energy cosmic rays, they must 
be photons. This is a problem because there is compelling evidence that the highest energy 
cosmic rays are not photons: 

1. The highest energy event observed by Fly's Eye is not likely to be a photon [|7|]. A 
photon of 300 EeV will interact with the magnetic field of the earth far above the atmosphere 
and disintegrate into lower energy cascades — roughly ten at this particular energy. The 
detector subsequently collects light produced by the fluorescence of atmospheric nitrogen 
along the path of the high-energy showers traversing the atmosphere. The anticipated 
shower profile of a 300 EeV photon is shown in Fig. 4. It disagrees with the data. The 
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FIG. 4: The composite atmospheric shower profile of a 3 x 10 eV gamma ray shower calculated 
with Landau-Pomeranchuk-Migdal (dashed) and Bethe-Heitler (solid) electromagnetic cross sec- 
tions. The central line shows the average shower profile and the upper and lower lines show 1 a 
deviations — not visible for the BH case, where lines overlap. The experimental shower profile is 
shown with the data points. It does not fit the profile of a photon shower. 

observed shower profile does fit that of a primary proton, or, possibly, that of a nucleus. 
The shower profile information is sufficient, however, to conclude that the event is unlikely 
to be of photon origin. 

2. The same conclusion is reached for the Yakutsk event that is characterized by a huge 
number of secondary muons, inconsistent with an electromagnetic cascade initiated by a 
gamma ray. 

3. The AGASA collaboration claims evidence for "point" sources above lOEeV. The 
arrival directions are however smeared out in a way consistent with primaries deflected by 
the galactic magnetic field. Again, this indicates charged primaries and excludes photons. 
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4. Finally, a recent reanalysis of the Haverah Park disfavors photon origin of the primaries 

Neutrino primaries are definitely ruled out. Standard model neutrino physics is under- 
stood, even for EeV energy The average x of the parton mediating the neutrino interaction 
is of order x ~ \Jm^Js ~ 10~ 6 so that the perturbative result for the neutrino-nucleus 
cross section is calculable from measured HERA structure functions. Even at 100 EeV a 
reliable value of the cross section can be obtained based on QCD-inspired extrapolations 
of the structure function. The neutrino cross section is known to better than an order of 
magnitude. It falls 5 orders of magnitude short of the strong cross sections required to make 
a neutrino interact in the upper atmosphere to create an air shower. 

Could EeV neutrinos be strongly interacting because of new physics? In theories with 
TeV-scale gravity, one can imagine that graviton exchange dominates all interactions and 
thus erases the difference between quarks and neutrinos at the energies under consideration. 
The actual models performing this feat require a fast turn-on of the cross section with energy 



that violates S-wave unitarity pi pi 15, 16, 17, 18, 19, 20, 21, 22, 23 



We have exhausted the possibilities. Neutrons, muons and other candidate primaries one 
may think of are unstable. EeV neutrons barely live long enough to reach us from sources 
at the edge of our galaxy. 



D. A Three Prong Assault on the Cosmic Ray Puzzle 

We conclude that, where the highest energy cosmic rays are concerned, both the ac- 
celerator mechanism and the particle physics are enigmatic. The mystery has inspired a 
worldwide effort to tackle the problem with novel experimentation in three complementary 
areas of research: air shower detection, atmospheric Cerenkov astronomy and underground 
neutrino astronomy. While some of the future instruments have additional missions, all are 
likely to have a major impact on cosmic ray physics. 



1. Giant Cosmic Ray Detectors 

With super-GZK fluxes of the order of a single event per square kilometer, per century, 
the outstanding problem is the lack of statistics; see Fig. 2 and Fig. 3. In the next five years, 
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a qualitative improvement can be expected from the operation of the HiRes fluorescence 
detector in Utah. With improved instrumentation yielding high quality data from 2 detectors 
operated in coincidence, the interplay between sky transparency and energy measurement 
can be studied in detail. We can safely anticipate that the existence of super-GZK cosmic 
rays will be conclusively demonstrated by using the instrument's calorimetric measurements. 
A mostly Japanese collaboration has proposed a next-generation fluorescence detector, the 
Telescope Array. 

The Auger air shower array is confronting the low rate problem with a huge collection 
area covering 3000 square kilometers on an elevated plain in Western Argentina. The in- 
strumentation consists of 1600 water Cerenkov detectors spaced by 1.5 km. For calibration, 
about 15 percent of the showers occurring at night will be viewed by 3 HiRes-style fluores- 
cence detectors. The detector is expected to observe several thousand events per year above 
10 EeV and tens above 100 EeV. Exact numbers will depend on the detailed shape of the 
observed spectrum which is, at present, a matter of speculation. 

2. Gamma rays from Cosmic Accelerators 

An alternative way to identify the source (s) of the highest energy cosmic rays is illustrated 
in Fig. 5. The cartoon draws our attention to the fact that cosmic accelerators are also cosmic 
beam dumps which produce secondary photon and neutrino beams. Accelerating particles to 
TeV energy and above requires relativistic, massive bulk flows. These are likely to originate 
from the exceptional gravitational forces associated with dense cores of exploding stars, 
inflows onto supermassive black holes at the centers of active galaxies, annihilating black 
holes or neutron stars. In such situations, accelerated particles are likely to pass through 
intense radiation fields or dense clouds of gas surrounding the black hole. This leads to 
the production of secondary photons and neutrinos that accompany the primary cosmic 
ray beam. An example of an electromagnetic beam dump is the UV radiation field that 
surrounds the central black hole of active galaxies. The target material, whether a gas of 
particles or of photons, is likely to be tenuous enough that the primary beam and the photon 
beam are only partially attenuated. However, shrouded sources from which only neutrinos 
can emerge, as in terrestrial beam dumps at CERN and Fermilab, are also a possibility. 

The astronomy event of the 21st century could be the simultaneous observation of TeV- 
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NEUTRINO BEAMS: HEAVEN & EARTH 
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FIG. 5: Diagram of cosmic accelerator and beam dump. See text for discussion. 

gamma rays, neutrinos and gravitational waves from cataclysmic events associated with the 
source of the highest energy cosmic rays. 

We first concentrate on the possibility of detecting high-energy photon beams. Af- 
ter two decades, ground-based gamma ray astronomy has become a mature science 



2% Pq , [26| , \27\ , |28| , [29| . A large mirror, viewed by an array of photomultipliers, collects 



the Cerenkov light emitted by air showers and images the showers in order to determine the 
arrival direction and the nature of the primary particle. These experiments have opened 
a new window in astronomy by extending the photon spectrum to 20 TeV, and possibly 
beyond. Observations have revealed spectacular TeV-emission from galactic supernova rem- 
nants and nearby quasars, some of which emit most of their energy in very short bursts of 
TeV-photons. 

But there is the dog that didn't bark. No evidence has emerged for the tt° origin of TeV 
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radiation. Therefore, no cosmic ray sources have yet been identified. Dedicated searches for 
photon beams from suspected cosmic ray sources, such as the supernova remnants IC433 
and 7-Cygni, came up empty handed. While not relevant to the topic covered by this paper, 
supernova remnants are theorized to be the sources of the bulk of the cosmic rays that are 
of galactic origin. However, the evidence is still circumstantial. 

The field of gamma ray astronomy is buzzing with activity to construct second-generation 
instruments. Space-based detectors are extending their reach from GeV to TeV energy with 
AMS and, especially, GLAST, while the ground-based Cerenkov collaborations are design- 
ing instruments with lower thresholds. Soon, both techniques should generate overlapping 
measurements in the 10 — 10 2 GeV energy range. All ground-based air Cerenkov experiments 
aim at lower threshold, better angular and energy resolution, and a longer duty cycle. One 
can, however, identify three pathways to reach these goals: 

1. larger mirror area, exploiting the parasitic use of solar collectors during nighttime 



(CELESTE, STACEY and SOLAR II) |0] 



2. better, or rather, ultimate imaging with the 17m MAGIC mirror, [EH 



3. larger field of view and better pointing and energy measurement using multiple tele- 
scopes (VERITAS, HEGRA and HESS). 



The Whipple telescope pioneered the atmospheric Cerenkov technique. VERITAS [32| is 



an array of 9 upgraded Whipple telescopes, each with a field of view of 6 degrees. These can 
be operated in coincidence for improved angular resolution, or be pointed at 9 different 6 
degree bins in the night sky, thus achieving a large field of view. The HEGRA collaboration 
|33| is already operating four telescopes in coincidence and is building an upgraded facility 
with excellent viewing and optimal location near the equator in Namibia. 

There is a dark horse in this race: Milagro ||34|| . The Milagro idea is to lower the threshold 



of conventional air shower arrays to 100 GeV by instrumenting a pond of five million gallons 
of ultra-pure water with photomultipliers. For time-varying signals, such as bursts, the 
threshold may be even lower. 
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3. Neutrinos from Cosmic Accelerators 



How many neutrinos are produced in association with the cosmic ray beam? The answer 
to this question, among many others [35. 3(J, provides the rational for building kilometer- 
scale neutrino detectors. 

Let's first consider the question for the accelerator beam producing neutrino beams at 
an accelerator laboratory. Here the target absorbs all parent protons as well as the muons, 
electrons and gamma rays (from tt° — > 7 + 7) produced. A pure neutrino beam exits the 
dump. If nature constructed such a "hidden source" in the heavens, conventional astronomy 
will not reveal it. It cannot be the source of the cosmic rays, however, for which the dump 
must be partially transparent to protons. 

In the other extreme, the accelerated proton interacts, thus producing the observed high- 
energy gamma rays, and subsequently escapes the dump. We refer to this as a transparent 
source. Particle physics directly relates the number of neutrinos to the number of observed 
cosmic rays and gamma rays ||37|| . Every observed cosmic ray interacts once, and only once, 
to produce a neutrino beam determined only by particle physics. The neutrino flux for such 



a transparent cosmic ray source is referred to as the Waxman-Bahcall flux pBL |39], HO, £11 



and is shown as the horizontal lines labeled "W&B" in Fig. 6. The calculations is valid for 
E ~ 100 PeV. If the flux is calculated at both lower and higher cosmic ray energies, however, 
larger values are found. This is shown as the non-flat line labeled "transparent" in Fig. 6. 
On the lower side, the neutrino flux is higher because it is normalized to a larger cosmic 
ray flux. On the higher side, there are more cosmic rays in the dump to produce neutrinos 
because the observed flux at Earth has been reduced by absorption on microwave photons, 
the GZK-effect. The increased values of the neutrino flux are also shown in Fig. 6. The 
gamma ray flux of tt° origin associated with a transparent source is qualitatively at the level 
of observed flux of non-thermal TeV gamma rays from individual sources [[37]]. 

Nothing prevents us, however, from imagining heavenly beam dumps with target densities 
somewhere between those of hidden and transparent sources. When increasing the target 
photon density, the proton beam is absorbed in the dump and the number of neutrino- 
producing protons is enhanced relative to those escaping the source as cosmic rays. For the 
extreme source of this type, the observed cosmic rays are all decay products of neutrons 
with larger mean-free paths in the dump. The flux for such a source is shown as the upper 
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log^/GeV) 

FIG. 6: The neutrino flux from compact astrophysical accelerators. Shown is the range of pos- 
sible neutrino fluxes associated with the the highest energy cosmic rays. The lower line, labeled 
"transparent", represents a source where each cosmic ray interacts only once before escaping the 
object. The upper line, labeled "obscured", represents an ideal neutrino source where all cosmic 
rays escape in the form of neutrons. Also shown is the ability of AMANDA and IceCube to test 
these models. 

horizontal line in Fig. 6. 

The above limits are derived from the fact that theorized neutrino sources do not overpro- 
duce cosmic rays. Similarly, observed gamma ray fluxes constrain potential neutrino sources 
because for every parent charged pion (7r ± — > I + u), a neutral pion and two gamma rays 
(7T° — > 7 + 7) are produced. The electromagnetic energy associated with the decay of neutral 
pions should not exceed observed astronomical fluxes. These calculations must take into ac- 
count cascading of the electromagnetic flux in the background photon and magnetic fields. 
A simple argument relating high-energy photons and neutrinos produced by secondary pions 
can still be derived by relating their total energy and allowing for a steeper photon flux as 
a result of cascading. Identifying the photon fluxes with those of non-thermal TeV photons 
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emitted by supernova remnants and blazers, we predict neutrino fluxes at the same level as 
the Waxman-Bahcall flux. It is important to realize however that there is no evidence that 
these are the decay products of 7r°'s. The sources of the cosmic rays have not been revealed 
by photon or proton astronomy 0, • 

For neutrino detectors to succeed they must be sensitive to the range of fluxes covered 
in Fig. 6. The AMANDA detector has already entered the region of sensitivity and is elim- 
inating specific models which predict the largest neutrino fluxes within the range of values 
allowed by general arguments. The IceCube detector, now under construction, is sensitive 
to the full range of beam dump models, whether generic as or modeled as active galaxies 
or gamma ray bursts. IceCube will reveal the sources of the cosmic rays or derive an upper 
limit that will qualitatively raise the bar for solving the cosmic ray puzzle. The situation 
could be nothing but desperate with the escape to top-down models being cut off by the 
accumulating evidence that the highest energy cosmic rays are not photons. In top-down 
models, decay products predominantly materialize as quarks and gluons that materialize as 
jets of neutrinos and photons and very few protons. We will return to top-down models at 
the end of this review. 

II. HIGH-ENERGY NEUTRINO TELESCOPES 
A. Observing High-energy Neutrinos 

Although details vary from experiment to experiment, high-energy neutrino telescopes 
consist of strings of photo-multiplier tubes (PMT) distributed throughout a natural 
Cerenkov medium such as water or ice. Typical spacing of PMT is 10-20 meters along 
a string with string spacing of 30-100 meters. Such experiments can observe neutrinos of 
different flavors over a wide range of energies by using a variety of methods: 

• Muon neutrinos that interact via charged current interactions produce a muon (along 
with a visible hadronic shower if the neutrino is of sufficient energy). The muon travels 
through the medium producing Cerenkov radiation which is detected by an array of 
PMT. The timing, amplitude (number of Cerenkov photons) and topology of the PMT 
signals is used to reconstruct the muon's path. The muon energy threshold for such a 
reconstruction is typically in the range of 10-100 GeV. 
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FIG. 7: Although IceCube detects neutrinos of any flavor, at TeV-EeV energies, it can identify 
their flavor and measure their energy in the ranges shown. Filled areas: particle identification, 
energy, and angle. Shaded areas: energy and angle. 

To be detected, a neutrino must interact via charged current and produce a muon 
with sufficient range to reach the detector. The probability of detection is therefore 
the product of the interaction probability (or the inverse interaction length A" 1 = na u ) 
and the range of the muon R^: 



P ~ 

1 V^>[1 — 



no v R u 



(11) 



where n is the number density of target nucleons, a v is the charged current interaction 



cross section [fj8| and the range is ~ 5 m per GeV for low energy muons. The muon 
range is determined by catastrophic energy loss (brehmsstrahlung, pair production 
and deep inelastic scattering) for muons with energies exceeding ~ 500 GeV |46| , 



Muon, tau or electron neutrinos which interact via charged or neutral current inter- 
actions produce showers which can be observed when the interaction occurs within 
or close to the detector volume. Even the highest energy showers penetrate water or 
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ice less than 10 m, a distance short compared to the typical spacing of the PMT. The 
Cerenkov light emitted by shower particles, therefore, represents a point source of light 
as viewed by the array. The radius over which PMT signals are produced is 250 m 
for a 1 PeV shower; this radius grows or decreases by approximately 50 m with every 
decade of shower energy. The threshold for showers is generally higher than for muons 
which limits neutral current identification for lower energy neutrinos. The probability 
for a neutrino to interact within the detector's effective area and to generate a shower 
within its volume is approximately given by: 

P u ^ ~ na u L (12) 

where a v is the charged+neutral current interaction cross section, L is the length of 
the detector along the path of the neutrino and n, again, is the number density of 
target nucleons. 

Tau neutrinos are more difficult to detect but produce spectacular signatures at PeV 
energies. The identification of charged current tau neutrino events is made by observ- 
ing one of two signatures: double bang events [[|9|, |50], ^TJ and lollypop events |5^, |53 



Double bang events occur when a tau lepton is produced along with a hadronic shower 
in a charged current interaction within the detector volume and the tau decays pro- 
ducing a electromagnetic or hadronic shower before exiting the detector (as shown in 
Fig. 8). Below a few PeV, the two showers cannot be distinguished. Lollypop events 
occur when only the second of the two showers of a double bang event occurs within 
the detector volume and a tau lepton track is identified entering the shower over several 
hundred meters. The incoming r can be clearly distinguished from a muon. A muon 
initiating a PeV shower would undergo observable catastrophic energylosses. Lollypop 
events are useful only at several PeV energies are above. Below this energy, tau tracks 
are not long enough to be identified. 

A feature unique to tau neutrinos is that they are not depleted in number by absorption 
in the earth. Tau neutrinos which interact producing a tau lepton generate another tau 
neutrino when the tau lepton decays, thus only degrading the energy of the neutrino 
H Hil |7|. 



Although MeV scale neutrinos are far below the energies required to identify individual 
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events, large fluxes of MeV electron anti-neutrinos interacting via charged current 
could be detected by observing higher counting rates of individual PMT over a time 
window of several seconds. The enhancement rate in a single PMT will be buried in 
dark noise of that PMT. However, summing the signals from all PMT over a short 
time window can reveal significant excesses, for instance form a galactic supernova. 

With these signatures, neutrino astronomy can study neutrinos from the MeV range to the 
highest known energies (~ 10 20 eV). 



B. Large Natural Cerenkov Detectors 

A new window in astronomy is upon us as high-energy neutrino telescopes see first 



light [SBfl. Although neutrino telescopes have multiple interdisciplinary science missions, 



the search for the sources of the highest-energy cosmic rays stands out because it most di- 



rectly identifies the size of the detector required to do the science [46], |47| . For guidance 
in estimating expected signals, one makes use of data covering the highest-energy cosmic 
rays in Fig. 2 and Fig. 3 as well as known sources of non-thermal, high-energy gamma rays. 
Estimates based on this information suggest that a kilometer-scale detector is needed to see 
neutrino signals as previously discussed. 

The same conclusion is reached using specific models. Assume, for instance, that gamma 
ray bursts (GRB) are the cosmic accelerators of the highest-energy cosmic rays. One can 
calculate from textbook particle physics how many neutrinos are produced when the particle 
beam coexists with the observed MeV energy photons in the original fireball. We thus 
predict the observation of 10-100 neutrinos of PeV energy per year in a detector with a 
square kilometer effective area. GRB are an example of a generic beam dump associated 
with the highest energy cosmic rays. We will work through this example in some detail in 
later sections. In general, the potential scientific payoff of doing neutrino astronomy arises 
from the great penetrating power of neutrinos, which allows them to emerge from dense 
inner regions of energetic sources. 

The strong scientific motivations for a large area, high-energy neutrino observatory lead to 
the formidable challenges of developing effective, reliable and affordable detector technology. 
Suggestions to use a large volume of deep ocean water for high-energy neutrino astronomy 
were made as early as the 1960s. Today, with the first observation of neutrinos in the Lake 
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FIG. 8: Simulation of an ultra high-energy tau lepton generated by the interaction of a 10 PeV tau 
neutrino (first shower), followed by the decay of the secondary tau lepton (second shower). The 
shading represents the time sequence of the hits. The size of the dots corresponds to the number 
of photons detected by the individual photomultipliers. 

Baikal and the South Pole neutrino telescopes, there is optimism that the technological 
challenges of building neutrino telescopes have been met. 

Launched by the bold decision of the DUMAND collaboration to construct such an in- 
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FIG. 9: The arrival times of the Cerenkov photons in 6 optical sensors determine the direction of 
the muon track. 

strument, the first generation of neutrino telescopes is designed to reach a large telescope 



area and detection volume for a neutrino threshold of order 10 GeV |59|, |6(J, |6T|]. This rela- 
tively low threshold permits calibration of the novel instrumentation on the known flux of 
atmospheric neutrinos. The architecture is optimized for reconstructing the Cerenkov light 
front radiated by an up-going, neutrino-induced muon. Up-going muons must be identified 
in a background of down-going, cosmic ray muons which are more than 10 5 times more 
frequent for a depth of ~l-2 kilometers. The earth is used as a filter to screen out the 
background of down-going cosmic ray muons. This makes neutrino detection possible over 
the hemisphere of sky faced by the bottom of the detector. 

The optical requirements on the detector medium are severe. A large absorption length is 
needed because it determines the required spacing of the optical sensors and, to a significant 
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extent, the cost of the detector. A long scattering length is needed to preserve the geometry 
of the Cerenkov pattern. Nature has been kind and offered ice and water as natural Cerenkov 
media. Their optical properties are, in fact, complementary. Water and ice have similar 
attenuation length, with the roles of scattering and absorption reversed. Optics seems, at 
present, to drive the evolution of ice and water detectors in predictable directions: towards 
very large telescope area in ice exploiting the long absorption length, and towards lower 
threshold and good muon track reconstruction in water exploiting the long scattering length. 



1. Baikal, ANT ARES, Nestor and NEMO: Northern Water 

Whereas the science is compelling, we now turn to the challenge of developing effective 
detector technology. With the termination of the pioneering DUMAND experiment, the 
efforts in water are, at present, spearheaded by the Baikal experiment |62], |63], 53], |65 



The Baikal Neutrino Telescope is deployed in Lake Baikal, Siberia, 3.6 km from shore at a 
depth of 1.1km. An umbrella-like frame holds 8 strings, each instrumented with 24 pairs 
of 37-cm diameter QUASAR photomultiplier tubes. Two PMT are required to trigger in 
coincidence in order to suppress the large background rates produced by natural radioactivity 
and bioluminescence in individual PMT. Operating with 144 optical modules (OM) since 
April 1997, the NT-200 detector was completed in April 1998 with 192 OM. Due to unstable 
electronics only ~ 60 channels took data during 1998. Nevertheless 35 neutrino-induced up- 
going muons were identified in the first 234 live days of data; see Fig. 10 for a 70 day sample. 
The neutrino events are isolated from the cosmic ray muon background by imposing a 
restriction on the chi-square of the fit of measured photon arrival times and amplitudes to a 
Cherenkov cone, and by requiring consistency between the reconstructed trajectory and the 
spatial locations of the OMs reporting signals. In order to guarantee a minimum lever arm 
for track fitting, they only consider events with a projection of the most distant channels on 
the track larger than 35 meters. This does, of course, result in a higher energy threshold. 
Agreement with the expected atmospheric neutrino flux of 31 events shows that the Baikal 
detector is understood. Stability and performance of the detector have improved in 1999 
and 2000 data taking |65 . 



The Baikal site is competitive with deep oceans, although the smaller absorption length 
of Cerenkov light in lake water requires a somewhat denser spacing of the OMs. This does, 
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FIG. 10: Angular distribution of muon tracks in the Lake Baikal experiment after the cuts described 
in the text. 

however, result in a lower threshold which is a definite advantage, for instance for oscillation 
measurements and WIMP searches. They have shown that their shallow depth of 1 kilometer 
does not represent a serious drawback. A significant advantage is that the site has a seasonal 
ice cover which allows reliable and inexpensive deployment and repair of detector elements. 

In the following years, NT-200 will be operated as a neutrino telescope with an effective 
area between 10 3 and 5 x 10 3 m 2 , depending on energy. Presumably too small to detect 
neutrinos from extraterrestrial sources, NT-200 will serve as the prototype for a larger 
telescope. For instance, with 2000 OMs, a threshold of 10 to 20 GeV and an effective 
area of 5 x 10 4 to 10 5 m 2 , an expanded Baikal telescope could fill the gap between present 
underground detectors and planned high threshold detectors of cubic kilometer size. Its key 
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advantage would be low energy threshold. 

The Baikal experiment represents a proof of concept for future deep ocean projects that 
have the advantage of larger depth and optically superior water. Their challenge is to 
find reliable and affordable solutions to a variety of technological challenges for deploying a 
deep underwater detector. Several groups are confronting the problem; both NESTOR and 
ANTARES are developing rather different detector concepts in the Mediterranean. 



The NESTOR collaboration |)6| |67|, |68fl , as part of a series of ongoing technology tests, 
is testing the umbrella structure which will hold the OMs. They have already deployed 
two aluminum "floors", 34m in diameter, to a depth of 2600m. Mechanical robustness 
was demonstrated by towing the structure, submerged below 2000 m, from shore to the site 
and back. These tests should soon be repeated with two fully instrumented floors. The 
cable connecting the instrument to the counting house on shore has been deployed. The 
final detector will consist of a tower of 12 six-legged floors vertically separated by 30 m. 
Each floor contains 14 OMs with four times the photocathode area of the commercial 8 inch 
photomultipliers used by AMANDA and ANTARES. 

The detector concept is patterned along the Baikal design. The symmetric up/down 
orientation of the OMs will result in uniform angular acceptance and the relatively close 
spacings will result in a low energy threshold. NESTOR does have the advantage of a 
superb site off the coast of Southern Greece, possibly the best in the Mediterranean. The 
detector can be deployed below 3.5 km relatively close to shore. With the attenuation length 
peaking at 55 m near 470 nm, the site is optically similar to that of the best deep water sites 
investigated for neutrino astronomy. 

The ANTARES collaboration |)9], [70], [HJ is currently constructing a neutrino telescope 
at a 2400 m deep Mediterranean site off Toulon, France. The site is a trade-off between 
acceptable optical properties of the water and easy access to ocean technology. Their detector 
concept requires remotely operated vehicles for making underwater connections. Results on 
water quality are very encouraging with an absorption length of 40 m at 467 nm and 20 m 
at 375 nm, and a scattering length exceeding 100 m at both wavelengths. Random noise, 
exceeding 50 khz per OM, is eliminated by requiring coincidences between neighboring OMs, 
as is done in the Lake Baikal design. Unlike other water experiments, they will point all 
photomultipliers sideways or down in order to avoid the effects of biofouling. The problem 
is significant at the Toulon site, but only affects the upper pole region of the OM. Relatively 
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weak intensity and long duration bioluminescence results in an acceptable deadtime of the 
detector. They have demonstrated their capability to deploy and retrieve a string, and have 
reconstructed down-going muons with 8 OMs deployed on the test string. 

The ANTARES detector will consist of 13 strings, each equipped with 30 stories and 
3 PMT per story. This detector will have an area of about 3 x 10 4 m 2 for 1 TeV muons 
- similar to AMANDA-II — and is planned to be fully deployed by the end of 2004. The 
electro-optical cable linking the underwater site to the shore was successfully deployed in 
October 2001. 

NEMO, a new R&D initiative based in Catania, Sicily has been mapping Mediterranean 
sites, studying mechanical structures and low power electronics. One hopes that with a suc- 
cessful pioneering neutrino detector of 10~ 3 km 3 in Lake Baikal and a forthcoming 10~ 2 km 3 
detector near Toulon, the Mediterranean effort will converge on a 10 _1 km 3 detector, possi- 



bly at the NESTOR site [|7^, [73||. For neutrino astronomy to become a viable science, several 
projects will have to succeed in addition to AMANDA. Astronomy, whether in the optical 
or in any other wave-band, thrives on a diversity of complementary instruments, not on "a 
single best instrument". 

2. AMANDA: Southern Ice 



Construction of the first-generation AMANDA-B10 detector |74|, [75|, [F7|, [78f was 



completed in the austral summer 96-97. It consists of 302 optical modules deployed at a 
depth of 1500-2000 m; see Fig. 11. Here the optical modules consist of 8-inch photomultiplier 
tubes and are controlled by passive electronics. Each is connected to the surface by a cable 
that transmits the high voltage as well as the anode current of a triggered photomultiplier. 
The instrumented volume and the effective telescope area of this instrument matches those 
of the ultimate DUMAND Octagon detector which, unfortunately, could not be completed. 

Depending on depth, the absorption length of blue and UV light in the ice varies between 
85 and 225 meters. The effective scattering length, which combines the mean-free path A 
with the average scattering angle 6 cts — (cosoyj ' "Metrics from 15 to 40 meters Because 



the absorption length of light in the ice is very long and the scattering length relatively 
short, many photons are delayed by scattering. In order to reconstruct the muon track, 
maximum likelihood methods are used, which take into account the scattering and absorption 
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FIG. 11: The AMANDA detector and a schematic diagram of an optical module. Each dot 
represents an optical module. The modules are separated by 20 meters in the inner strings 1-4, 
and by 10 meters in the outer strings 5-10. 



of photons as determined from calibration measurements [73]. A Bayesian formulation of 
the likelihood |80| , which accounts for the much larger rate of down-going cosmic-ray muon 
tracks relative to up-going signal, has been particularly effective in decreasing the chance 
for a down-going muon to be misreconstructed as up-going. 

Other types of events that might appear to be up-going muons must also be considered 



29 



and eliminated. Rare cases, such as muons which undergo catastrophic energy loss, for 
instance through bremsstrahlung, or that are coincident with other muons, must be inves- 
tigated. To this end, a series of requirements or quality criteria, based on the characteristic 
time and spatial pattern of photons associated with a muon track and the response of the 
detector, are applied to all events that, in the first analysis, appear to be up-going muons. 
For example, an event which has a large number of optical modules hit by photons unscat- 
tered (relative to the expected Cerenkov times of the reconstructed track) has a high quality. 
By making these requirements (or "cuts" ) increasingly selective, they eliminate more of the 
background of false up-going events while still retaining a significant fraction of the true up- 
going muons, i.e., the neutrino signal. Two different and independent analyses of the same 
data covering 138 days of observation in 1997 have been undertaken. These analyses yielded 
comparable numbers of up-going muons (153 in analysis A, 188 in analysis B). Comparison 
of these results with their respective Monte Carlo simulations shows that they are consistent 
with each other in terms of the numbers of events, the number of events in common, and, 
as discussed below, the expected properties of atmospheric neutrinos. 

In Fig. 12, from analysis A, the experimental events are compared to simulations of back- 
ground and signal as a function of the (identical) quality requirements placed on the three 
types of events: experimental data, simulated up-going muons from atmospheric neutrinos, 
and a simulated background of down-going cosmic ray muons. For simplicity in presentation, 
the levels of the individual types of cuts have been combined into a single parameter repre- 
senting the overall event quality, and the comparison is made in the form of ratios. Fig. 12 
shows events for which the quality level is 4 and higher. As the quality level is increased 
further, the ratios of simulated background to experimental data and experimental data to 
simulated signal both continue their rapid decrease, the former toward zero and the latter 
toward unity. Over the same range, the ratio of experimental data to the simulated sum of 
background and signal remains near unity. At an event quality of 6.9 there are 153 events 
in the sample of experimental data and the ratio to predicted signal is 0.7. The conclu- 
sions are that (1) the quality requirements have reduced the events from misreconstructed 
down-going muons in the experimental data to a negligible fraction of the signal and that 
(2) the experimental data behave in the same way as the simulated atmospheric neutrino 
signal for events that pass the stringent cuts. They estimate that the remaining signal is 
contaminated by instrumental background at 15 ± 7 percent. 
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FIG. 12: Reconstructed muon events in AMANDA-BIO are compared to simulations of background 
cosmic ray muons (BG MC) and simulations of atmospheric neutrinos (Signal MC atm v) as a 
function of "event quality", a variable indicating the severity of the cuts designed to enhance the 
signal. Note that the comparison is made in the form of ratios. 

The estimated uncertainty on the number of events predicted by the signal Monte Carlo 
simulation (which includes uncertainties in the high-energy atmospheric neutrino flux, the 
sensitivity of the optical modules, and the precise optical properties of the ice) is +40% 
to —50%. The observed ratio of experiment to simulation (0.7) and the expectation (1.0) 
therefore agree within errors. 

The shape of the zenith angle distribution from analysis B is compared to a simulation 
of the atmospheric neutrino signal in Fig. 13 in which the two distributions have been 
normalized to each other. The variation of the measured rate with zenith angle is reproduced 
by simulation to within the statistical uncertainty. Note that the tall geometry of the 
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FIG. 13: Reconstructed zenith angle distribution for AMANDA-BIO. The points mark the data 
and the shaded boxes a simulation of atmospheric neutrino events. The widths of the boxes indicate 
the error bars. The overall normalization of the simulation has been adjusted to match the data. 

detector strongly influences the dependence on zenith angle in favor of more vertical muons. 

Estimates of the energies of the up-going muons (based on simulations of the number of 
optical modules that participate in an event) indicate that the energies of these muons are 
in the range from 100 GeV to ~ 1 TeV. This is consistent with their atmospheric neutrino 
origin. 

The agreement between simulation and experiment shown in Fig. 12 and 13, taken to- 
gether with other comparisons of measured and simulated events, leads us to conclude that 
the up-going muon events observed by AMANDA are produced mainly by atmospheric neu- 
trinos. 

The arrival directions of the neutrinos observed in both analyses are shown in Fig. 14. 
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FIG. 14: Distribution in declination and right ascention of the up-going AMANDA-BIO events on 
the sky. 

A statistical analysis indicates no evidence for point sources in this sample. An estimate of 
the energies of the up-going muons indicates that all events have energies consistent with 
an atmospheric neutrino origin. This corresponds to a level of sensitivity to a diffuse flux 
of high-energy extra-terrestrial neutrinos of order dN/dE u = 10~ 6 E~ 2 cm -2 s _1 sr _1 GeV -1 , 
assuming an E~ 2 spectrum [^TJ. This upper limit excludes a variety of theoretical models 
which assume the hadronic origin of TeV photons from active galaxies and blazars. Searches 
for neutrinos from gamma ray bursts, magnetic monopoles, and for a cold dark matter signal 
from the center of the Earth yield limits comparable to or better than those from smaller 
underground neutrino detectors that have operated for a much longer period. 

Data are being taken now with the larger array, AMANDA-II consisting of an additional 
480 OMs. 



3. IceCube: A Kilometer-Scale Neutrino Observatory 

The IceCube project Q at the South Pole is a logical extension of the research and 
development work performed over the past several years by the AMANDA Collaboration. 
The optimized design for IceCube is an array of 4800 photomultiplier tubes each enclosed in a 
transparent pressure sphere to comprise an optical module similar to those in AMANDA. In 
the IceCube design, 80 strings are regularly spaced by 125 m over an area of approximately 
one square kilometer, with OMs at depths from 1.4 to 2.4 km below the surface. Each 
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string consists of OMs connected electrically and mechanically to a long cable which brings 
OM signals to the surface. The array is deployed one string at a time. For each string, a 
enhanced hot- water drill melts a hole in the ice to a depth of about 2.4 km in less than 2 
days. The drill is then removed from the hole and a string with 60 OMs vertically spaced 
by 17 m is deployed before the water re- freezes. The signal cables from all the strings are 
brought to a central location which houses the data acquisition electronics, other electronics, 
and computing equipment. 

Each OM contains a 10 inch PMT that detects individual photons of Cerenkov light 
generated in the optically clear ice by muons and electrons moving with velocities near the 
speed of light. 

Background events are mainly down-going muons from cosmic ray interactions in the 
atmosphere above the detector. The background is monitored for calibration purposes and 
background rejection by the IceTop air shower array covering the detector. 

Signals from the optical modules are digitized and transmitted to the surface such that a 
photon's time of arrival at an OM can be determined to within less than 5 nanoseconds. The 
electronics at the surface determines when an event has occurred (e.g., that a muon traversed 
or passed near the array) and records the information for subsequent event reconstruction 
and analysis. 

At the South Pole site (see Fig. 15), a computer system accepts the data from the event 
trigger via the data acquisition system. The event rate, which is dominated by down-going 
cosmic ray muons, is estimated to be 1-2 kHz. The technology that will be employed in 
IceCube has been developed, tested, and demonstrated in AMANDA deployments, in labo- 
ratory testing, and in simulations validated by AMANDA data. This includes the instrument 
architecture, technology, deployment, calibration, and scientific utilization of the proposed 
detector. There have been yearly improvements in the AMANDA system, especially in 
the OMs, and in the overall quality of the information obtained from the detector. In the 
1999/2000 season, a string was deployed with optical modules containing readout electronics 
inside the OM. The information is sent digitally to the surface over twisted-pair electrical 
cable. This option eliminates the need for optical fiber cables and simplifies calibration of 
the detector elements. This digital technology is the baseline technology of IceCube. For 
more details, see Ref. 



The construction of neutrino telescopes is overwhelmingly motivated by their discovery 
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FIG. 15: The South Pole site, showing the residential dome and associated buildings, the ski- 
way where planes land, the dark sector with the Martin A. Pomerantz Observatory in which the 
AMANDA electronics are housed, and a rough outline of where IceCube strings are to be placed. 

potential in astronomy, astrophysics, cosmology and particle physics. To maximize this 
potential, one must design an instrument with the largest possible effective telescope area 
to overcome the neutrino's small cross section with matter, and the best possible angular 
and energy resolution to address the wide diversity of possible signals. 

At this point in time, several of the new instruments (such as the partially deployed Auger 
array, HiRes, Magic, Milagro and AMANDA II) are less than one year from delivering results. 
With rapidly growing observational capabilities, one can realistically hope, almost 100 years 
after their discovery, the puzzling origin of the cosmic rays will be deciphered. The solution 
will almost certainly reveal unexpected astrophysics or particle physics. 

C. EeV Neutrino Astronomy 

At extremely high energies, new techniques can be used to detect astrophysical neutrinos. 
These include the detection of acoustic and radio signals induced by super-EeV neutrinos 
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interacting in water, ice or salt domes, or the detection of horizontal air showers by large 
conventional cosmic ray experiments such as the Auger array. 

Horizontal air showers are likely to be initiated by a neutrino because showers induced by 
primary cosmic rays are unlikely to penetrate the ~ 36,000g/cm 2 of atmosphere along the 
horizon. Isolated penetrating muons may survive but they can be experimentally separated 
from a shower initiated by a neutrino close to the detector. Horizontal air shower experiments 
can also use nearby mountains as a target, e.g. to observe the decay of tau leptons produced 
in charged current interactions in the moutain. The sensitivity of an air shower array 
to detect an ultra high-energy neutrino is described by its acceptance, expressed in units 
of km 3 water equivalent steradians (km 3 wesr). Typically only showers with zenith angle 
greater than ~ 70 degrees can be identified as neutrinos. This corresponds to a slant depth 
of ~ 2000 g/cm 2 . 

The acceptance of present air shower experiments, such as AGASA, is ~ 1km 3 we sr above 
10 10 GeV, and significantly less at lower energies. Auger will achieve ten times greater accep- 
tance at 10 9 GeV and 50 times greater near 10 12 GeV. Nitrogen fluorescence experiments 
also have the capability to detect neutrinos as nearly horizontal air showers with space- 
based experiments such as EUSO and OWL extending the reach of Auger. At this point 
we should point out however that the actual event rates of these experiments are similar 
to those for IceCube. Although IceCubes energy resolution saturates at EeV energies, the 
neutrinos are still detected with rates competitive with the most ambitious horizontal air 
shower experiments; for a more detailed comparison see Ref. p5, ^ . 



Radio Cerenkov experiments detect the Giga-Hertz pulse radiated by shower electrons 
produced in the interaction of neutrinos in ice. Also, the moon, viewed by ground-based 
radio telescopes, has been used as a target |37] . Above a threshold of ~ 1 PeV, the large 
number of low energy(~ MeV) photons in a shower will produce an excess of electrons over 
positrons by removing electrons from atoms by Compton scattering. These are the sources 
of coherent radiation at radio frequencies, i.e. above ~ 100 MHz. The mechanism is now 
well understood. The characteristics and the power of the pulses have been measured by 



dumping a photon beam in sand ||88|| . The results agree with calculations |p9[ . 

While many proposals exist, the most extensive effort to develop a radio neutrino detector 
is RICE (Radio Ice Cerenkov Experiment), which is located in the shallow ice above the 



AMANDA detector [90]. It consists of an 18-channel array of radio receivers distributed 
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within a 8 x 10 6 m 3 volume. The receivers, buried in the ice at depths of 100-300 meters, 
are sensitive over the range of 0.2-1 GHz, roughly corresponding to electron neutrinos with 
energy of several PeV and above. The ANITA collaboration proposes to fly a balloon-borne 
array of radio antennas on a circular flight over Antarctica. ANITA will detect earth- 
skimming neutrinos producing signals emerging from the ice along the horizon P2"| . 
With higher threshold but also greater effective area than RICE (about 1 million km 2 ), 
ANITA should be sensitive to GZK neutrinos after a lucky 30 day flight (or 3 normal flights 
of 10 days). 

EeV neutrino-induced showers can also be detected by acoustic emission resulting from 
local heating of a dense medium. Existing arrays of hydrophones, built in the earth's oceans 
for military application, could be used for the hydro-acoustic detection of neutrinos with 
extremely high energies; for a recent review see |)3] . 



III. COSMIC NEUTRINO SOURCES 



A List of Cosmic Neutrino Sources 



We have previously discussed generic cosmic ray producing beam dumps and their asso- 
ciated neutrino fluxes. We now turn to specific sources of high-energy neutrinos. The list of 
proposed sources is long and includes, but is not limited to: 

• Gamma Ray Bursts (GRB) 

GRB, outshining the entire universe for the duration of the burst, are perhaps the 



best motivated source for high-energy neutrinos [53, [95|, [96|]. Although we do not yet 
understand the internal mechanisms that generate GRB, the relativistic fireball model 
provides us with a successful phenomenology accommodating observations. It is very 
likely that GRB are generated in some type of cataclysmic process involving dying 
massive stars. GRB may prove to be an excellent source of neutrinos with energies 
from MeV to EeV and above. As we shall demonstrate further on, their fluxes can be 
calculated in a relatively model independent fashion. 

Other Sources Associated with Stellar Objects 

Other theorized neutrino sources associated with compact objects include supernova 
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remnants exploding into the interstellar medium |47], ^8], [59|], X-ray binaries |46|, 
TOOj TO ggg , microquasars 0, [10|, |ig and even the sun g6|, g7|, |105], fD§|, any of 



which could provide observable fluxes of high-energy neutrinos. 
Active Galactic Nuclei (AGN): Blazars 

Blazars, the brightest objects in the universe and the sources of TeV-energy gamma 
rays, have been extensively studied as potential neutrino sources. Blazar flares with 
durations ranging from months to less than an hour, are believed to be produced by 
relativistic jets projected from an extremely massive accreting black hole. Blazars may 
be the sources of the highest energy cosmic rays and, in association, provide observable 
fluxes of neutrinos from TeV to EeV energies. 

Neutrinos Associated with the Propagation of Cosmic Rays 

Very high-energy cosmic rays generate neutrinos in interactions with the cosmic mi- 
crowave background [ 107 , |108|| . This cosmogenic flux is among the most likely sources 
of high-energy neutrinos, and the most straightforward to predict. Furthermore, cos- 
mic rays interact with the Earth's atmosphere ||109| , |110|| and with the hydrogen con- 
centrated in the galactic plane [H5|, |111| , |112| , |113|| producing high-energy neutrinos. 
It has also been proposed that cosmic neutrinos themselves may produce cosmic rays 
and neutrinos in interactions with relic neutrinos u+Ub — > Z. This is called the Z-burst 



mechanism [114, 115, 116, 117, 118 



Dark Matter, Primordial Black Holes, Topological Defects and Top-Down Models 

The vast majority of matter in the universe is dark with its particle nature not yet 
revealed. The lightest supersymmetric particle, or other Weakly Interacting Massive 
Particles (WIMPs) propsed as particle candidates for cold dark matter, should become 
gravitationally trapped in the sun, earth or galactic center. There, they annihilate 
generating high-energy neutrinos observable in neutrino telescopes [|119| , |120| , |121| , |122| , 
123| , |124| , |125| . Another class of dark matter candidates are superheavy particles with 



GUT-scale masses that may generate the ultra high-energy cosmic rays by decay or 
annihilation, as well as solve the dark matter problem. These will also generate a 

127] , |128j , |13U[ |. Extremely high-energy neutrinos are 



substantial neutrino flux 



also predicted in a wide variety of top-down scenarios invoked to produce cosmic 
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rays, including decaying monopoles, vibrating cosmic strings ||131|, |132|| and Hawking 



radiation from primordial black holes ||133| , |134| , |135|| 



Any of these sources may or may not provide observable fluxes of neutrinos. History 
testifies to the fact that we have not been particularly successful at predicting the 
phenomena invariably revealed by new ways of viewing the heavens. We do, however, 
know that cosmic rays exist and that nature accelerates particles to super-EeV energy. 
In this review we concentrate on neutrino fluxes associated with the highest energy 
cosmic rays. Even here the anticipated flux depends on our speculation regarding the 
source. We will work through three much-researched examples: GRB, AGN and decays 
of particles or defects associated with the GUT-scale. The myriad of speculations 



have been recently reviewed by Learned and Mannheim ||47 |. We concentrate here on 



neutrino sources associated with the observed cosmic rays and gamma rays. 



B. Gamma Ray Bursts: A Detailed Example of a Generic Beam Dump 

1. GRB Characteristics 

Although there is no such thing as a typical gamma ray burst, observations of GRB 
indicate the following common characteristics: 

• GRB are extremely luminous events, often releasing energy of order one solar mass in 
gamma rays. Typically, L 7 ~ 10 51 to 10 54 erg/s is released over durations of seconds 
or tens of seconds. GRB are the most luminous sources in the universe. 

• GRB produce a broken power-law spectrum of gamma rays with 7 oc E~ 2 for E 1 > 
0.1-1 MeV and 7 oc E' 1 for £ 7 < 0.1-1 MeV |3§, pTj] , 



GRB are cosmological events. Redshifts exceeding z=4 have been measured [138|, |139 



GRB are rare. During it's operation, BATSE observed on average 1 burst per day 
within its field of view (~ 1/3 of the sky). Assuming that the rate of GRB does not 
significantly change with cosmological time, this corresponds to one burst per galaxy 
per million years. If GRB are beamed, they may be more common. 
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GRB produce afterglows of less energetic photons which extend long after the initial 
burst Pd, PH|, P2|. 



The durations of GRB follow a bimodal distribution with peaks near two seconds and 
20 seconds, although some GRB have durations ranging from milliseconds to 1000 
seconds | 143 |. Variations in the spectra occur on the scale of milliseconds [|143| , 144 
shown in Fig. 16 ||145|| . GRB afterglows can extend for days [ |1 13|| . 
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FIG. 16: An example of the temporal structure of a GRB as recorded by BATSE [145]. Note the 
two time scales: a duration of several seconds and a fluctuation time scale of a fraction of a second. 



2. A Brief History of Gamma Ray Bursts 

Gamma Ray Bursts (GRB) were accidentally discovered in the late 1960's by the military 
Vela satellites, intended to monitor nuclear tests in space forbidden by the Outer Space 
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Treaty between the United States and the Soviet Union ||146|| . The Vela observation of a 



short, intense burst of MeV gamma rays was originally considered to be a possible signal 
from an advanced extra-terrestrial civilization . The idea was quickly reconsidered. In 1973, 
the discovery was announced to the public [ |146[ | . Shortly after, the observation was verified 



by the Soviet IMP-6 satellite flUT 



Until the 1990 's, the high intensity of GRB led astronomers to the belief that they were 
galactic in origin. In 1991, the BATSE (Burst And Transient Satellite Experiment) detector 
on the Compton Gamma Ray Observatory was launched. BATSE observed roughly one 
burst per day within its field of view of about one third of the sky. The observations 



showed total isotropy of GRB over the entire sky, thus ruling out galactic origin ||148| . The 
cosmological origin of GRB implies that they release up to a solar mass of energy, in seconds 
time. Their cosmological origin was subsequently confirmed by afterglow observations, first 



made in 1997 by the Beppo-SAX satellite ||140|| . Afterglow observations were made in X-ray, 



optical and longer wavelengths with an angular resolution of arc-minute precision and with 
measurement of the redshift. To date, dozens of GRB afterglows have been observed, nearly 



all of which have resulted in the identification of the host galaxy [ 149 , 15C. 151 . 

Although progress has been made in our understanding of GRB, many questions remain 
unanswered. Most importantly, the progenitor(s) of GRB remain an open question. In the 
next section, we describe some of the most likely candidates. 

3. GRB Progenitors? 

The observed characteristics of GRB require an original event with a large amount of 
energy (~ M Q ) in a very compact volume (i?o ~ 100 km). The phenomenology that describes 
observations is that of a fireball expanding with highly relativistic velocity, powered by 
radiation pressure. The nature of the "inner engine" that initiates the fireball remains an 
open question. Afterglow observations have recently shown that GRB are predominantly 
generated in host galaxies and are likely the result of a stellar process. Research into a 
variety of stellar progenitors has been pursued. 

The "Collapsar" scenario, where a super-massive star undergoes core collapse resulting 
in a failed supernova, is one of the most common models proposed for the fireball's inner 



engine ||152| , |153|| . As matter falls into the black hole created in this process, gravitational 
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energy is transferee! to bulk kinetic energy and the fireball is generated. 

The strength of magnetic fields and the angular momentum of the stellar object (s) in- 
volved can play an important role in the dynamics of the core collapse process. For example, 
"Magnetars" are a subset of the core collapse model which result in a rapidly spinning neu- 
tron star with an extremely strong magnetic field [133, |155| , |156| . Objects with sufficient 
angular momentum can undergo a "supranova" process where their core collapse takes place 
in two stages, possibly separated by months or years [157 , 158|| . In this scenario, the ob- 
ject's large angular momentum prevents a fraction of the matter from falling into the fireball 
initially. 

Compact objects in close binary orbits are also likely candidates for fireball progenitors. 
The "hypernovae" scenario is similar to the core collapse models, but includes a secondary 
stellar object in the dynamics [159, 160 , 161 , 162j| . Similarly, neutron star binaries or neutron 
star-black hole binaries (or possibly white dwarf - neutron star or black hole binaries) which 
lose sufficient angular momentum through gravitation radiation can undergo a merger. Such 
a merger is expected to generate a black hole surrounded by debris. As this debris is accreted 



into the black hole, the required fireball is generated ||161|, |163|, |164|, |165|, |166|, |167|, |168 



Finally, if primordial strange hadrons exist, a "seed" of strange matter may start a 
chain reaction converting a neutron star into a strange star made entirely of strange matter 
169| , |1 70 , |171[ |. This conversion would release the majority of the star's binding energy as 



it contracts, thus generating a compact fireball similar to that required for GRB dynamics. 

Recent evidence indicates the presence of emission lines in GRB ||172|| . This evidence 
strengthens the argument for progenitors involving collapsing stars. 

The problem of GRB progenitors is likely to have an experimental solution. Possible 
progenitor-specific signatures may be found by using gravitational waves [|173|| or neutrinos 



as astronomical probes, or by more detailed study of afterglows [174, 175 



It is interesting to note that the bimodal distribution of GRB durations may be an 
indication of multiple GRB classes and associated progenitors. 

4- Fireball Dynamics 



a. The Fireball The dynamics of a gamma ray burst fireball is similar to the physics 
of the early universe. Initially, there is a radiation dominated soup of leptons and photons 
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and few baryons. This perfect fluid has the equation of state P = p/3 and is initially hot 
enough to freely produce electron-positron pairs. The luminosity of a burst can be related 
to the number density of photons n 7 : 

L = AnRlcn^E^, (13) 

where Ro is the initial radius of the source, i.e. prior to expansion. The optical depth of a 
photon before pair production is determined by the photon density and the interaction cross 
section |176| : 



R n La Th L 7 \ /100km\ ZlMeV\ 

Here Ai n t is the interaction length of a photon as a result of pair production and Thomson 
scattering. These cross sections are roughly equal with the Thomson cross section o"Th — 
10- 24 cm 2 . 

With an optical depth of order ~ 10 15 , photons are trapped in the fireball. This results in 



the highly relativistic expansion of the fireball powered by radiation pressure ||168| , |177|| . The 
fireball will expand with increasing velocity until it becomes transparent and the radiation is 
released. This results in the visual display of the GRB. By this time, the expansion velocity 
has reached highly relativistic values of order 7 ~ 300. 

Besides leptons and photons, the fireball contains some baryons. During expansion, 
the opaque fireball cannot radiate and any nucleons present are accelerated as radiation 
is converted into bulk kinetic energy. When the radiation is emitted, there is a transition 
from radiation to matter dominance of the fireball. At this stage, the radiation pressure is 
no longer important and the expanding fireball coasts without acceleration. The expansion 
velocity remains constant with 7 ~ 77 = that is determined by the amount of baryonic 
matter present, often referred to as the baryon loading | |1 78" , |1 79| | . The phenomenology will 



reveal values of 77 between 10 2 and 10 3 [|136| , |176| , 180]. The formidable appearance of the 



GRB display is simply associated with the large boost between the fireball and the observer 
who detects highly boosted energies and contracted times. 

The exploding fireballs original size, i?o, is that of the compact progenitor, for instance 
the black hole created by the collapse of a massive star. As the fireball expands the flow is 
shocked in ways familiar from the emission of jets by the black holes at the centers of active 
galaxies or mini-quasars. (A way to visualize the formation of shocks is to imagine that 
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Fireball Frame Observer Frame 




FIG. 17: Diagram of GRB fireball kinematics assuming no beaming. Primed quantities refer to 
the comoving frame. Unprimed quantities refer to the observer's frame. 

infalling material accumulates and chokes the black hole. At this point a blob of plasma 
is ejected. Between these ejections the emission is reduced.) The net result is that the 
expanding fireball is made up of multiple shocks. These are the sites of the acceleration of 
particles to high-energy and the seeds for the complex millisecond structures observed in 
individual bursts; see Fig. 16. (Note that these shocks expand with a range of velocities and 
they will therefore collide providing a mechanism to accelerate particles to high-energy.) The 
characteristic width of these shocks in the fireball frame is SR' = 'ycAt, where At ~ 0.01 sec. 



For an alternative scenario; see Ref. [181, 182, 183 



An expanding shock is seen by the observer as an expanding shell of thickness cAt = R$ 
and radius R; see Fig. 17. Here At is the time scale of fluctuations in the burst fireball; it is 
related to R' by: 

R! = 7 2 cAt = 7 2 i?o, (15) 
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with primed quantities referring from now on to the frame where the fireball is at rest. Two, 
rather than a single 7-factor, relate the two quantities because of the geometry that relates 
the radius R' to the time difference between photons emitted from a shell expanding with a 
velocity v; see Fig. 17 [fL37||. Introducing the separation d of the two photons along the line 



of sight, we note that 



d n ,«, R' ^ v 2 . R' 

At # _ ~ 1 _ -) ~ — , (16) 

c c c Ic c A cy z 



using the relativistic approximation that 1 + ^ ~ 1. 

We next calculate the energy of the burst. In the observer frame 

E = U xV oc T 4 x R 2 AR oc T 4 R 2 oc 7 4 T /4 J R 2 , (17) 

where {7 is the energy density and V is the fireball shell volume. In the fireball frame, 

E' = U' x V oc T' 4 R' 3 oc 7 6 T' 4 . (18) 

Energy conservation requires that E and E' remain constant during expansion of the fireball. 
In the fireball frame, this results in the usual blackbody relation that T' 4 is proportional to 
R'~ 3 or, using Eq. 18, proportional to 7 -6 . Substituting into the expression for E, we obtain 
that 

E oc -f~ 2 R 2 , (19) 

or, because E is constant, that 

i?oc7. (20) 

Thus we obtain the important result that, with expansion, the 7-factor grows linearly with 
R until reaching the maximum value 77. 

b. The Observed GRB Spectrum: Synchroton and Inverse Compton Scattering The 
broken power-law gamma ray spectrum of GRB, with two distinct spectral slopes, is far from 
a blackbody spectrum. The observed spectrum, therefore, clearly indicates that fireball pho- 
tons do not sufficiently interact to thermalize prior to escaping the fireball. After escaping, 
the photons show spectral features characteristic of the high-energy, non-thermal emission by 
supernova remnants and active galaxies. Here photons up to MeV-energy can be produced 
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by synchrotron radiation, with some reaching, possibly, up to TeV energies by inverse Comp- 
ton scattering on accelerated electrons |34], |136|| . These processes have been modeled for the 



expanding fireball and successfully accommodate observed GRB spectra. This represents a 
major success of the relativistic fireball phenomenology |184j , |185| , [TSR 187, 188|. 



To produce the non-thermal spectrum, special conditions must prevail [ |136| , |137| , |176| , |180 
The photons must not thermalize prior to the time when the shock becomes transparent and 
the observed radiation released. Conversely, if they decouple too early, there is insufficient 
time for the synchotron and inverse Compton scattering processes to produce the observed 
spectrum. This requires that the expansion time of the shockwave and the time Ai nt /c for 
photons and electrons to interact by Thompson scattering be similar: 

— ~ Kcctt)- 1 . (21) 

7c 

Here ot is the Thompson cross section and n e is the electron number density in the fireball. 
The latter can be related to the mass flux M with the assumption that n e = n p : 

M = 4irR 2 cyn p m p . (22) 

As required by mass conservation, M is independent of R since n p oc 7 -3 and R oc 7. In 
terms of luminosity, 

ML 
He ^ — (23) 

4:nR 2 c~fm p ^i]c^ATiR 2 m p 
where rj = L/Mc 2 is the ratio of luminosity to mass previously introduced, rj is also referred 
to as the dimensionless entropy and, as previously derived, 7 ~ 77 after expansion of the 
fireball. We can rewrite the condition of Eq. 21 for producing the observed non-thermal 
spectrum as 

Rn e a T La T 1 77 * 4 



7 7 2 77c 3 47ri?m p 7 s rj 
where the critical dimensionless entropy 77* is defined as 



, Lo T \ 1/4 / L \ ^/lOO km\ l / A 

"=U^) *™*{wwi) (— ) • (25) 

At decoupling, rj = 7 and Eq. 24 is satisfied, provided 77 = 77*. The condition for the fireball to 
produce the correct non-thermal, synchrotron/inverse Compton spectrum is realized with the 
expansion time matching the Thompson scattering time. For values of 77 that are significantly 
larger (smaller), the decoupling of the radiation will occur too early (late) thus limiting rj , 
as well as the final value of 7, to the range of 10 2 — 10 3 . 
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(24) 



c. Jets and Beaming Observations imply that the total amount of energy emitted in 
gamma rays by a GRB are typically in the range of 10 52 — 10 54 ergs, i.e. a large fraction of 
a solar mass. For some bursts, it may exceed a solar mass. This, as well as the difficulty 
of converting such an unusually large fraction of primary energy into gamma rays, strongly 
suggests that GRB are beamed. Beaming reduces the total amount of energy by a factor 
of Q/4tt, where Q is the solid angle in which the observed gamma rays are emitted. Most 
proposed progenitors naturally predict a rotating stellar source that is likely to produce 
beamed emission. Relativistic beaming is possible down to an angular size of Q > 7~ 2 , 
although larger angles are of course possible [|189| , |190| , |191| . 



In the presence of beaming, the number of bursts is increased by a factor of 4ir/Q in order 
to account for bursts that do not point towards earth and are, therefore, not observed. For 
typical Lorentz factors of 7 ~ 300 and a minimum beaming angle of Q > 7~ 2 ~ 10~ 5 , on the 



order of one GRB per galaxy per year is required to accommodate the observations ||190 



It is important to note that most of the diffuse neutrino fluxes calculated in this review are 
independent of beaming because the reduced energy for a single burst is compensated by 
their increased frequency. 



5. Ultra High-energy Protons From GRB? 



As previously discussed, it may be possible to accelerate protons to energies above 10 20 eV 
in GRB shocks pO], [192], [L93|. GRB within the GZK radius of 50-100 Mpc, could therefore 
be the source of the ultra high-energy cosmic rays (UHECR's) |96|, |97|, |180| , |192| , |193| , |19 1| . 



195| , [L96| , |197| , |198|| . To accelerate protons to this energy, several conditions have to be 
satisfied. First, the acceleration time t a ~ AR^/c, where A is a factor of order 1 and 
Rl = E/eBj is the Larmor radius, must not exceed the duration of the burst R/^c, 



AE R 

■ • 7 ' 



'jeB 



(26) 



or 



AE , / E \/10 n m 
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(27) 
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Second, energy losses due to synchrotron radiation must not exceed the energy gained by 
acceleration. The synchrotron loss time is given by 

tsyn = = • (28) 

c cn e or 

The number density of electrons (in the rest frame) is given by 



Therefore, 



67T Quince 3 
tsyn = o T ml&Bi lv = a T mlEB*' { } 



For synchrotron energy losses to be less than the energy gained by acceleration, 

6nmtc 3 E , x 

^^x— -, (31) 



(TtitiIEB 2 ~ jecB' 
or 

1 / ry \ 2 /10 20 eV\ 2 

B 4 >< (^r) ■ < 32 > 

Combining above requirements, we get 

-Ki^)(^)-- B 4---(^ 2 (^) 2 ' w 

or 

fl^xlO'Veters^)'^) 3 . (34) 
From simple fireball kinematics, we previously derived that 

R < l 2 cAt, (35) 
where At is ~ 10 msec. Combining this with Eq. 34 leads to the final requirement: 

- 42xlo '°(T) 2 (lo^v) 3sflS72cAt ' (36) 



or 



48 



which can indeed be satisfied for the values of 7 = 10 2 — 10 3 previously derived from fireball 
phenomenology. We conclude that bursts with Lorentz factors > 100 can accelerate protons 
to ~ 10 20 eV. The long acceleration time of 10-100 seconds implies however that the fireball 
extends to a large radius where surrounding matter may play an important role in the 
kinematics of the expanding shell. 

Finally, it can be shown that proton energy losses from p — 7 interactions will not interfere 
with acceleration to high-energy. These will, in fact, be the source of high-energy neutrinos 
associated with the beam of high-energy protons. We will discuss this further on. 

6. Neutrino Production in GRB: the Many Opportunities 

Several mechanisms have been proposed for the production of neutrinos in GRB. We 
summarize them first: 

• Thermal Neutrinos: MeV Neutrinos 

As with supernovae, GRB are expected to radiate the vast majority of their initial 
energy as thermal neutrinos. Although the details are complex and are likely to 
depend on the progenitor, a neutrino spectrum with a higher temperature than a 
supernova may be expected. Observation is difficult because of the great distances to 
GRB, although we should keep in mind that a nearby GRB may not be less frequent 
than a galactic supernova; see the section on beaming. 

• Shocked Protons: TeV-EeV Neutrinos 

Protons accelerated in GRB can interact with fireball gamma rays and produce pions 
that decay into neutrinos. While astronomical observations provide information on the 
fireball gamma rays, the proton flux is a matter of speculation. A definite neutrino 
flux is however predicted when assuming that GRB produce the highest energy cosmic 
rays. 

• Decoupled Neutrons: GeV Neutrinos 

In a GRB fireball, neutrons can decouple from protons in the expanding fireball. If 
their relative velocity is sufficiently high, their interactions will be the source of pions 
and, therefore, neutrinos. Typical energies of the neutrinos produced are much lower 
than those resulting from interactions with gamma rays. 
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7. Thermal Me V Neutrinos from GRB 



As is the case for a supernova, we expect that in GRB, thermal neutrinos are produced 



escaping with the majority of the total energy ||199| , |200|| . The dynamics are somewhat 
different, however. The temperature of the photons in a supernova is of order lOMeV as 
derived from the familiar estimate: 



U 52 er K s/ I R ) ' 1 j 



V 4htJ U0 52 ergs 

where cr ~ 5.67 x 10 _8 km/sec 3 is the Stefan-Boltzman constant and Uj is the energy density 
of the initial plasma. We copy this estimate for the neutrinos produced in the pre-fireball 
phase of a GRB taking into account the much larger energy emitted, typically 2 orders of 
magnitude: 

Here h u = 2x3x7/8 = 21/4 and hj = 2 are the degrees of freedom available to each 
particle type. This yields neutrinos with average energy: 

Ew 8 V^/iookmX 3 / 4 



3.15T,, ~ 90MeV 



10 54 ergs 

The flux of neutrinos can now be calculated from the average energy of roughly 100 Me V 
per individual neutrino and the total energy available: 

k ^ g x 10l0km _ 2 ( _^_^V_R_N'/Y3000Mpcy 

E VtWe 4nD 2 Vl0 54 ergs7 V 100 km J \ D J v ' 

Neutrinos with this energy are below threshold for the detection methods previously de- 
scribed. For supernova 1987A, predominantly electron anti-neutrinos were observed by the 
electromagnetic showers generated inside the detector by the process z7 e +p — > n + e + . Under- 
ground detectors are too small to detect the above flux because of the cosmological distance 
to the source. In large under-ocean detectors, the signal is drowned in the ~ 50 kHz noise 
in the photomultipliers from potassium decay. Only a large Cerenkov detector embedded in 
sterile ice can possibly detect GRB but even here the signal-to-noise is marginal unless the 



source is within our local cluster. Detailed calculations have been performed in Ref. [199 
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8. Shocked Protons: Pe V Neutrinos 



Assuming that GRB are the sources of the highest energy cosmic rays and that the 
efficiency for conversion of fireball energy into the kinetic energy of protons is similar to 
that for electrons, the production of PeV neutrinos is a robust prediction of the relativistic 
fireball model [S3, 137| , 194 1 . Neutrinos are produced in interactions of accelerated protons 



with fireball photons, predominantly via the processes 

P7 -> A -> nvr ± (42) 

and 

P7 -> A -> p7T° (43) 

which have very large cross sections of 10 _28 cm 2 . The charged 7r's subsequently decay pro- 
ducing charged leptons and neutrinos, while neutral 7r's may generate high-energy photons 
observable in TeV energy air Cerenkov detectors. 

For the center-of-mass energy of a proton-photon interaction to exceed the threshold 
energy for producing the A-resonance, the comoving proton energy must exceed 

m\ — ml 



4b; 



Therefore, in the observer's frame, 



*>"^(&)'(T)- 

resulting in a neutrino energy 

E„4 W£( >7 X 10»,V(±)'(^). (46) 

Here (x p ^ n ) ~ .2 is the average fraction of energy transferred from the initial proton to the 
produced pion. The factor of 1/4 is based on the estimate that the 4 final state leptons in 
the decay chain ^ — > u^n — > v^ei> e v e equally share the pion energy 

We already discussed the fireballs expansion and the formation of shocks that are the 
sites of the acceleration of particles to high-energy and the seeds for the complex millisecond 
structures observed in individual bursts. The characteristic width of a shock in the fireball 
frame is AR' = 7cAt, where At ~ 0.01 sec. 
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As the kinetic energy in fireball protons increases with expansion, a fraction of this energy 
is converted into pions once the protons are accelerated above threshold for pion production. 
The fraction of energy converted to pions is estimated from the number of proton interactions 
occurring within a shock of characteristic size AR': 

AR' 

U - -T—(Xp-*n)- (47) 

The proton interaction length A P7 in the photon fireball is given by 

— = n 7 (TA. (48) 

Ap 7 

Here n 7 is the number density of photons in the fireball frame and a a ~ 10 _28 cm 2 is the 
proton-photon cross section at the A-resonance. The photon number density is the ratio of 
the photon energy density and the photon energy in the comoving frame: 

7 E> 7 WrvARi)/ \ 7 J' 1 yj 

Using the fireball kinematics of Eqs. 15 and 16 

Thus we obtain the fraction of proton energy converted to 7r's in the expansion: 

f ~ L 7 1 ^a(x^) ^ / L, \ / lMeV \ /SOOxV-OlsecN 

U ~ E^^At 4ttc 2 " vl0 52 erg/sA E 1 A 7 A A* )' 1 ] 

For L ~ 10 52 erg/sec, At ~ 10 msec and 7 ~ 300, this fraction is on the order of 10 
percent. This quantity strongly depends on the Lorentz factor 7. Even modest burst-to- 
burst fluctuations in 7 around the average value of 300 can result in a PeV neutrino flux 
dominated by a few bright bursts; we will return to a discussion of fluctuations further on 
20H g02| . 



In order to normalize the neutrino flux we introduce the assumption that GRB are the 
source of cosmic rays above the ankle of the cosmic ray spectrum near ~ 3 x 10 18 eV 



180, 192, 193, 194, 195, M 197, 198]. The flux in neutrinos can then be simply 



obtained from the total energy injected into cosmic rays and the average energy of a single 



neutrino 94 



4tt EL ~ 4vr E v ~ 4vr E v \2 U h dt ) ' [b ] 
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or 



0, = 2 x 10-W VV- f 7Xl ° eV ) ( -U f * ) f f/ dt 3 J) , (53) 

r V E„ / V0.125/ VlOGyr/ Ux 10 44 Mpc'V" 1 / 

where f# ~ 10 Gyrs is the Hubble time and dE/dt ~ 4 x 10 44 ergs Mpc~ 3 yr -1 is the injection 

rate of energy into the universe in the form of cosmic rays above the ankle. 

PeV neutrinos are detected by observing a charged lepton produced in the charged current 

interaction of a neutrino near the detector. For instance, the probability to detect a muon 

neutrino within a Cerenkov neutrino telescope's effective area is given by Eq. 11. At TeV- 

PeV energies the function P u ^^ can be approximated by 

P v ^ ~ 1.7 x lO-X-l(TeV), (54) 

where E Uj0bs = E v /{1 + z) is the observed neutrino energy. The rate of detected events is 
the convolution of the flux with the probability of detecting the neutrino 

£max Jjji 

N e vents= I " (f) u P u ^-^- ^ 25 knT 2 yr -1 (55) 

J ^thresh 

This rate is significantly enhanced when fluctuations in distance, energy and (possibly) 
Lorentz factor are considered. The event rate is likely to be dominated by a few bright 
bursts rather than by a diffuse flux. 

With the ability to look for GRB neutrino events in coincidence with gamma ray obser- 
vations, i.e. in short time windows over which very little background accumulates, there is 
effectively no background for this neutrino signature of GRB. 



9. Stellar Core Collapse: Early Te V Neutrinos 

The core collapse of massive stars is, arguably, the most promising mechanism for gener- 
ating GRB. The fireball produced is likely to be beamed in jets along the collapsed object's 
rotation axis. The mechanism is familiar from observations of jets associated with the cen- 
tral black hole in active galaxies. The jets subsequently run into the stellar matter accreting 
onto the black hole. If the jets successfully emerge from the stellar envelope a GRB results. 
Interestingly, failed "invisible" jets which do not emerge will not produce a GRB display 



but will still produce observable neutrinos |[203| , |205 | 



A beamed GRB jet expanding with a Lorentz factor 7j et ~ 100 — 1000 through the 
stellar envelope, will be slowed down resulting in a smaller Lorentz factor at its leading edge 
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7f ^ 7jct- Therefore, the fast particles in the tail will catch up with the slow particles in the 
leading edge and collide with a Lorentz factor 7 7j Ct /27f by simple addition of relativistic 
velocities. Once the jet emerges from the infalling stellar matter around R ~ 10 6 km, the 
density drops to around ~ 10~ 7 g/cm 3 . Matching the energy densities on either side of the 
shock front requires 

7f 2 x l(T 7 g/cm 3 ~ (7jet/27f) n p m p , (56) 

where n p is the comoving proton number density in the jet. It is related to the luminosity 
of the burst; see Eq. 23 : 

Up = 4irR^ t m pC f (57) 

Here L iso = is the inferred from a non-beamed flux. Eqs. 56 and 57 determine the value 
of 7f: 

_ / L iso y/ 4 _ 3 x / L iso y/Y io-ycm^ vY io'km y/ 2 



7f "li6^^J " 3x Ito^J { — - P — ) {-it) ■ (58) 

With this low value of the Lorentz factor, the fireball remains opaque to gamma rays as 
described in the section on fireball dynamics. The radiation thermalizes with a temperature 
determined by its energy density U' and the Stefan-Boltzman law U' = AaT" 1 / c. We find 

7 v 4a J vio- 7 g/cm 3 y v 3 y y ' 

The rest of the calculation follows the previous section. Protons traveling through this 
thermal photon plasma produce pions, predominantly via the A-resonance, for energies 

B ;>^, 7 x^ G ev(^X). (60, 

Primed energies refer to the comoving frame with 7f£" = E. Neutrinos emerge from the 
interactions with energy 

^4<^)^10Tev(f)(^|X). (61) 

For mildly relativisitic conditions, the fraction of protons converted to pions is expected 
to be high in the very dense plasma, i.e. of order unity. The neutrino flux observed from a 
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single GRB at a distance D is calculated from the total energy emitted in neutrinos and the 
average energy of a single neutrino: 

E iso (x p ^) , _ 2 ( E iso V lOTeV y 3000Mpc \ 2 

This flux of TeV neutrinos from a single burst results in 

Agents ~ 4>u X ~ -05 — (63) 



. 10 53 ergs/ V -D 

events observed per year, per square kilometer of the detector. Here E vobs = E v j (1 + z) is 
the observed neutrino energy and, for TeV-energy neutrinos, we used the approximation 

P v ^ ~ 1.3 x 10- 6 ^, o6s (TeV). (64) 

The event rate is low. Bursts within a few hundred megaparsecs (~ 10 bursts per year as 
well as an additional unknown number of "invisible" bursts from failed GRB) may produce 
multiple TeV neutrino events in a kilometer scale detector. This signature is unique to 
supernova progenitors. 

10. UHE Protons From GRB: EeV Neutrinos 

Recent observations of GRB afterglows show evidence that GRB explode into an inter- 
stellar medium, consistent with the speculations that they are collapsing or merging stars. 
Shocks will be produced when the GRB runs into the interstellar medium, including a re- 
verse shock that propagates back into the burst ejecta. Electrons and positrons in the 
reverse shock radiate an afterglow of eV-keV photons that represent a target for neutrino 
production by ultra high-energy protons accelerated in the burst |204) . 



The fraction of proton energy going into 7r-production is calculated as before following 
Eq. 47, 

fn ^ -T— fcp-w), (65) 

U ~ E : -y 4 At Aire 2 ' [ } 

J - , -y,min Irs *±/iu 

where 7 rs is the Lorentz factor of the reverse shock. For the afterglow, the relevant time 
scale is 10-100 seconds and the luminosity is L 7 oc E~ x l 2 |176[| . E ljmin , the minimum photon 
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energy to produce pions via the A-resonance, is given by: 

Therefore 10 21 eV protons can kinematically produce 7r's on photons with energy as low as 
10 eV. Combining Eq. 66 and Eq. 67 we find 

/10eV\/ E v \V2 /300\ 5 /20sec\ , . 

Note that above keV energy, the photon luminosity follows the broken spectrum with L 1 oc 
E~ l and, therefore, f n oc E p rather than f w oc E^ 2 . 

Associating the accelerated beam with the observed ultra high-energy cosmic ray flux, 
dN p /dE p ~ AE~ 2 , where A ~ 5 x lO^m-W^GeV, and using E v ~ .05£ p , see Eq.46, 
the resulting neutrino flux is given by 
dN dN 

ii m ~ dE {Ep = 20Eu) x fAEp = 20Ev) ' (69) 

dN V/ „. A . „ , 2 / lOeV \ / 20^ \V2 /300\ 5 /20sec\ 

— («~Ax(20«-x.003 x( )(— ) (-)(- Sr ). (70) 



^ 2 * ^ x ( ™ ) (^(^-VV^v. (71 ) 

It is important to note that if the burst occurs in a region of higher density gas, as can 
be the case for a collapsing star, reverse shocks are produced earlier and, therefore, with 
smaller Lorentz factors. This results in /„. ~ 1. Then, 

^■(Eu) * ^(E P = 20E U ) xl~ix (20£,)- 2 , (72) 

^E 2 U (E U ) ~ 10~ 6 m- 2 s- 1 sr- 1 GeV. (73) 

In either case, the result is only valid above the threshold energy required to generates pions 
via the A-resonance, 

Below this threshold, ultra high-energy protons may still interact with non-thermal MeV 
photons, however. 
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The event rate in a neutrino telescope is calculated following Eq. 11. In the high-energy 
appoximation, 

P v ^ ~ 1.2 x 10- 2 £°;l(EeV). (75) 

This yields 

/•5xl0 10 GeV 

N events ~ 10~ 6 E- 2 x 3.7 x lO-^^d^nrV^sr- 1 ~ .01km- 2 yr -1 - (76) 

J7xl0 8 GeV 

This is a very small rate indeed. The neutrino energy is, however, above the threshold for 
EeV telescopes using acoustic, radio or horizontal air shower detection techniques. This 
mechanism may represent an opportunity for detectors with very high threshold, but also 
large effective area to do GRB physics. 



11. The Decoupling of Neutrons: GeV Neutrinos 

The conversion of radiation into kinetic energy in the fireball will accelerate neutrons along 
with protons, especially if the progenitor involves neutron stars. Protons and neutrons are 
initially coupled by nuclear elastic scattering. If the expansion of the fireball is sufficiently 
rapid the neutrons and protons will no longer interact. Neutrons decouple from the fireball 
while protons are still accelerated. Protons and neutrinos may then achieve relative velocities 
sufficient to generate pions which decay into GeV neutrinos | |206| , |207| | . We define the ratio 
of neutrons to protons as, 

Z = ~, (77) 

n p 

which initially remains constant during expansion. The fraction of neutrons which generate 
pions is calculated in the same way as in Eq. 47, 

/. - (78) 
We can relate the density of nucleons to the density of photons by the dimensionless entropy, 



r]m p c z 

Following the arguments used in our discussion of PeV neutrinos, we arrive at 



7^7 



(79) 



f hi ^£ (go) 

Jn m p ?n 3 A£47rc 2 (l + 0' 
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where a np ~ 3 x 10 cm is the neutron-proton cross section for pion production. As the 
neutrons and protons decouple, f w approaches unity. Using the fact that 7 asymptotically 
approaches rj at the end of expansion, we see that decoupling occurs for 

„>„ =( ^ y /4 ~ 4 Qor — y/yiookmy/y^y^ (81) 

In fact, the requirement for exceeding the threshold for 7r production is 77 > 1.2r) np [2"0~6 



The scattering time is therefore longer than the expansion time by a factor 1.2 4 ~ 2.1 and 
1 — e -2 ' 1 = 88% of the neutrons scatter. If threshold were exceeded by, say rj = 1.5r) np , then 
more than 99% of the neutrons would scatter. It is therefore a reasonable approximation 
to assume that all neutrons produce 7r's as long as r\ is above threshold. The number of 
neutrons in the fireball is large with 

m p c 2 l + £r) V10 53 ergs/ U + f/ V rj J K J 

Above the pion threshold, every neutron interacts with a proton producing one of the fol- 
lowing: 

• p + n — > p + p + n~ — > + fi~ — > e~ + v e + + ^ 

• p + n — > n + ra + 7r + — > + /i + — ► e + + z/ e + + f M 

• p + n— >p + n, + 7r ^7-|-7 

Thus on average, each interaction produces two 30-50 MeV muon neutrinos and two 30-50 
MeV electron neutrinos or two 70 MeV photons. The observed neutrino energy is 

£ vo6s ~ 30 - 50MeVx ~ 6 - 10 GeV x (-L)(—L-\ (83) 

i + z Uoo/u + z/ v ; 

This energy is below the threshold of neutrino telescopes with the possible exception of 
Baikal and ANTARES provided it is built with a sufficiently dense arrangement of the 
photomultipliers. 

P V ^{E V ~ GeV) ~ 10- 7 E,, obs (GeV) (84) 

parametrizes the chance of detecting a ~GeV muon neutrino in ANTARES as described by 
Eq. 11. This leads to an event rate from a single burst of 
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events — N v ^^.j~)2 ^ u ^^ (®^) 

where N v = 2N n is the number of muon neutrinos emitted by the burst, A e fj is the detector's 
effective area and D is the distance to the burst. Positioning, for simplicity, all bursts at 
z=l, this reduces to 



N events ~ N bursts x ^2 iV n ^^P^ M (8GeV)) (86) 

VlOOOyr-V V10 53 ergs/ Vl + C/ V2 + 2-2 v / rT^/ V.lkm 2 / J V ; 

This estimate is optimistic and the rate quite small. Observation requires a burst with 
favorable fluctuations in distance and energy. We discuss this important aspect of GRB 
detection next. 



12. Burst-To-Burst Fluctuations and Neutrino Event Rates 

We have focused so far on the diffuse flux of neutrinos produced collectively by all GRB. 
Given the possibility to observe individual GRB, it is important to consider the fluctuations 
that may occur in the dynamics from burst to burst. In the presence of large fluctuations, 
the relevant observable becomes the number of neutrinos observed from individual bursts 



201| , |202|| . We will consider fluctuations in the distance to the burst D, and in its energy 
E. One can further contemplate fluctuations in the Lorentz factor, 7, or equivalently in 77, 
the dimensionless entropy, although it is not clear whether energy and Lorentz factor are 
independent quantities. 

The number of neutrinos from a single source at a distance D is proportional to D~ 2 . 
Considering burst-to-burst fluctuations in distance will enhance the fluxes previously calcu- 
lated. This can be see as follows 

/ rfH — dN{D)D- 2 _ / ctH — dD _ Dl ve ct KuhUe Dl 



avc 



j*****dN(D)D£ tf**** B*D£dD |(ct Hu bbic) 3 >t Hu bbie) 2 ' (88) 
where we have assumed an isotropic Euclidian distribution of sources. In contrast, the 

average distance of a burst D ave is 

/ D 2 dD = - / D 2 dD -> 7J avc = 2- 1 / 3 ct H ubbie. (89) 

JO I Jo 
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Therefore, the event rate is enhanced by spatial fluctuations by a factor 



3 x 2~ 2/3 ~ 1.9. 



(90) 



This enhancement factor should be applied to the diffuse flux. It also represents an in- 
creased probability to observe a single burst yielding multiple neutrino events. Assuming a 
more realistic cosmological distribution for starburst galaxies increases the effect of spatial 
fluctuations by an additional factor of 2 above the Euclidian result. 

Next, we consider energy. Observations show that approximately one out of ten bursts is 
ten times more energetic than the median burst and approximately one out of a hundred is 
a hundred times more energetic. This results in an enhancement of the neutrino signal by a 
factor of 



If instead of a step function, a smooth distribution j£L oc E~ 2 covering two orders of mag- 
nitude is used, the enhancement factor is increased to 4.6. 

We thus reach the very important conclusion that correct averaging of bursts over their 
cosmological spatial distribution and over their observed energy distribution enhances the 
neutrino signal by approximately one order of magnitude. Variations in distance and energy 
affect all GRB neutrino fluxes previously discussed, regardless of the production mechanism. 

The third, and most important, quantity to vary from burst-to-burst is the Lorentz 
factor, or entropy. Its variation can greatly modify some fluxes, e.g. as ~ 7 4 for PeV 
neutrinos or ~ 7 5 for EeV neutrinos, and it barely affects others, e.g. GeV neutrinos 
produced by decoupled neutrons. The range over which 7 can vary is, however, limited 
|136| , \137[ |176| , |180|| . The Lorentz factor i) must be roughly 300 to generate the observed 
non-thermal gamma ray spectrum, ii) must be large enough to produce the highest energy 
cosmic rays of ~ 10 20 eV, and iii) the energy of shocked protons must be sufficient to be 
above threshold for producing pions on fireball gamma rays. 

The last condition requires that 



.89 X Median + -10 X 10 £ med ian + -01 X 100 £ med ian 



2.9 E, 



median • 



(91) 




(92) 



and the second that 




(93) 
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following Eq. 37. Once again, A is a factor of order unity, E the maximum proton energy 
generated and At the time scale associated with the rapid variations observed in single 
bursts. 

Together, these constraints limit any variation of the Lorentz factor to at most one order 
of magnitude. This does, however, correspond to a variation over four (five) orders of 
magnitude in the fraction of energy converted to 7r's which yield PeV (EeV) neutrinos. A 
change in Lorentz factor modifies the peak neutrino energy of the neutrinos as 7 2 and is, 
therefore, a secondary effect relative to the fraction of energy converted to pions. 

How the Lorentz factor of GRB vary from burst-to-burst is an open question. It is inter- 
esting to note that even with conservative fluctuations of a factor of two or so (say between 
200 and 400), the difference in neutrino flux exceeds one order of magnitude. Combined with 
variations in distance, and perhaps burst energy, the occasional nearby and bright (small 7, 
high E, or both) burst becomes a superior experimental signature to the diffuse flux. More 
detailed modeling has been presented in Ref. ||201] , [202j| . For an alternative perspective, see 



reference ||208 



13. The Effect of Neutrino Oscillations 

Perhaps the most important discovery of particle physics in the last decade is the oscil- 
lation of neutrinos. The fact that neutrinos can change flavor as they propagate can have 
an effect on the neutrino fluxes observed on earth. For propagation over cosmological dis- 
tances, the neutrino survival probablities become simple. For the example of neutrinos from 
pion decay, where there is an initial ratio of 1 to 2 to for electron, muon and tau flavors, 
respectively, we expect approximately a ratio of 1 to 1 to 1 after oscillations. Therefore 
most of the event rates calculated for muon neutrinos have been overestimated by a factor 
of two. At sufficient energies, however, electron and tau neutrinos can be observed as well 
thus counteracting this effect in part. 

For a review of neutrino oscillations, see Ref. [ |209| |. 
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C. Blazars: the Sources of the Highest Energy Gamma rays 



1. Blazar Characteristics 



Active Galactic Nuclei (AGN) are the brightest sources in the universe. They are of 
special interest here because some emit most of their luminosity at GeV energy and above. 
A subset, called blazers, emit high-energy radiation in collimated jets pointing at the earth. 
They have the following characteristics: 

• Although less luminous than GRB, with inferred isotropic luminosities of ~ 10 45 — 
10 49 ergs, they radiate this luminosity over much longer time periods with regular flares 
extending for days pT|, gll], HI. The energetics require a black hole roughly one 
billion times more massive than our sun. 

• Blazars produce radiation from radio waves to TeV gamma rays with enhancements 
in E 2 dN/dE or uF(v) in the IR to X-ray and the MeV-TeV range |2E| |14], H ED1 > 
217|| . Roughly 60 sources have been observed in the MeV-GeV range by the EGRET 
instrument on the Compton Gamma ray Observatory. A handful of TeV observations 
have been reported thus far |210|, |2T|, [ZTg |220|, ^ |222|. EGRET has found that 



blazers produce a typical 7 oc E spectrum in the MeV-GeV range. This spectrum 
may extend above or below this range ||210| , ETTl |223| ]. There appears to be an inverse 
relationship between blazar luminosity and peak emission. High luminosity blazars 
tend to peak in the GeV and optical bands | |211| , |215| , |224| , |225| , while low luminosity 
blazars tend to peak in the TeV and X-ray bands |m|, |2|, |2|, |27|, |228| . 



Blazars are cosmological sources. They appear less distant than GRB only because 
their lower luminosity makes distant sources difficult to observe ||210 , 



The time scale of flaring in blazar luminosity varies from a fraction of a day to years. 
This range of time scales indicates a sub-parsec engine; cAt ~ R' — > R' ~ 10~ 3 — 10 _1 



pc 210, 211, 223, 229 . 



Prior to the launch of EGRET, the only survey above 100 MeV energy had been made 
by the COS-B satellite. The COS-B survey revealed the first extragalactic gamma ray 
source, the active galaxy 3C-272 ||230|| . After EGRET was launched in 1991, several gamma 
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ray blazars were discovered. The majority of these were flat-spectrum radio-loud quasars, 
although some were BL Lac objects which emit into the TeV range ||223| , |229| , [231|| . 



Three of the EGRET sources have also been observed by ground-based atmospheric 
Cerenkov telescopes previously discussed. Future observations of blazars will include lower- 
threshold (below 50 GeV) ground-based telescopes, as well as next generation satellites such 
as GLAST |2T{], gg, g3g. 

2. Blazar Models 

It is widely believed that blazars are powered by accreting supermassive black holes with 
masses of ~ 10 7 M© or more. Some of the infalling matter is reemitted and accelerated in 
highly beamed jets aligned with the rotation axis of the black hole. 

It is generally agreed upon that synchrotron radiation by accelerated electrons is the 
source of the observed IR to X-ray peak in the spectrum ||211| , |231| , |233| , |234j , |235|| . Inverse 



Compton scattering of synchrotron or, possibly, other ambient photons by the same electrons 
can accommodate all observations of the MeV-GeV second peak in the spectrum. There is a 
competing explanation for the second peak, however. In hadronic models R211| , [236| , |237| , |238| , 



239], |240| , [241] , |242|| , MeV-GeV gamma rays are generated by accelerated protons interacting 



with gas or radiation surrounding the black hole. Pions produced in these interactions decay 
into the observed gamma rays and not-yet observed neutrinos. This process is accompanied 
by synchrotron radiation of the protons. 

The basic dynamics of the blazar is common to all models. Relativistic jets are generated 
with substructure that takes the form of "blobs" or "sheets" of matter traveling along the 
jet with Lorentz factors of 10-100. As previously discussed, in order to accommodate the 
observation of flares, the thickness of these sheets must be less than jcAt ~ 10~ 2 parsec, 
much smaller than their width, which is of the order of 1 parsec. It is in these blobs that 
shocks produce TeV gamma rays and high-energy neutrinos. 

Several calculations of the neutrino flux from active galactic nuclei have been performed 
137j |243| , g4|, HH B B B HI HI Hill IHH HI Hi- In the following sections, we 



describe two examples that illustrate the mechanisms for neutrino production in hadronic 
blazars. 
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FIG. 18: Diagram of Blazar Kinematics 
3. Highly Shocked Protons: Ee V Blazar Neutrinos 

If protons are present in blazar jets, they may interact with photons via the A-resonance 
to generate pions which then decay into very high-energy neutrinos. This process is similar 
to the process generating PeV and EeV neutrinos in GRB. There are some important dif- 
ferences, however. First, the Lorentz factor of the motion of the blob, traveling towards the 
observer, is typically smaller than for GRB shells. It can be constrained by considering the 
energy carried by the highest energy gamma rays observed in blazars. These gamma rays 
escape the blob and must, therefore, be below the energy threshold for pair production with 
ambient photons whose energy typically peaks around ~ 10 eV; the UV bump. Contem- 
plating the observations of gamma rays above 15 TeV in Markarian 501 |[218| , |219|| , evading 
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pair production requires 

# 7 ,max# 7 ,peak < T^e; ( 94 ) 
/ p \ 1/2 / p , \ 1/2 

7> 25(^^1 (f^E) . (95) 

1 U5TeW VIOeVy 1 ; 

We will therefore consider Lorentz factors in the range of 10 to 100. 

A second difference between blazars and GRB is the geometry of the shocked material. 

Instead of a shell, blobs can be treated as roughly spherical. The energy density is 

U' = ^= L ^ At = ^ (96) 

7 InR' 3 7 |7r( 7 cAt) 3 47rc 3 7 4 At 2 { ' 



Except for geometry, this is identical to Eq. 49. We obtain a number density of photons 

% 3L„ 

Following the arguments leading to Eqs. 51, we obtain a rather large conversion of energy 
into pions 



H ~< E> iE.nc^At 2 (97) 



K_ „ L 7 1 3a A (x p ^) ^ / L 7 X ZlOeVX /30X V lOOOsec X 

U ~ A P7 " E, 7 2 At 4ttc 2 ~ ^ x { m ^ eTg/s ) { E , \ At J 

When f n approaches unity, pions will be absorbed before decaying into neutrinos, requiring 
the substitution of f n by 1 — e - ^. Note that fluctuations in 7 are not as important as 
for GRB because the fraction of energy transferred to pions varies as 7 2 rather than 7 4 . 
Moreover, it quickly saturates near unity, especially for relatively low 7, i.e. bright flares or 
short time scales. 

For protons to photoproduce pions on photons with the ubiquitous UV photons of ~ 10 eV 
energy, 



m a — ml 



Therefore, in the observer's frame, 



If blazars are the sources of the highest energy cosmic rays, protons are accelerated to this 
energy and will generate accompanying neutrinos with energy 

E v = \(x p ^)E p > 7 x 10 17 ev(^) 2 (^). (101) 
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The neutrino flux from blazars can be calculated in the same way as for GRB 

^iM^-^iy--^ < io2) 

„ = w-^-v. (L^£L\ (*) (^L) {1 _ e -, )e a-.-, (103) 



E„ J UOGyr/ U x 10 44 / 

using Eq. 52 and Eq. 53. This a flux on the order of 100 km _2 yr _1 over 2tt steridian. The 

number of detected events is obtained from Eqs. 11, 

2 i P x 10 17 eVy/2 / t H \( dE/dt \ . . 

Events ~ foPu^u. ~ 10 knT 2 yr — =— tt 1 g 7 104 

y ^ V Ej, / VlOGyr/ V4 x 10 44 Mpc- 3 yr^V V ; 

for f w ~ .35. For values of 7 varying from 10 to 100, the number of events varies from 1 
to 70 events km _2 yr _1 , respectively. Observation in a kilometer-scale detector should be 
possible. 

The greatest uncertainty in this calculation is associated with the requirement that blazars 
accelerate protons to the highest observed energies. In the next section we present an 
alternative mechanism for producing neutrinos in blazars which does not invoke protons of 
such high-energy. 



4- Moderately Shocked Protons: TeV Blazar Neutrinos 

In line-emitting blazers, external photons with energies in the keV-MeV energy range are 
known to exist. They are clustered in clouds of quasi-isotropic radiation. Protons of lower 
energy relative to those contemplated in the previous section, can photoproduce pions in 



interactions with these clouds ||250| , [251| , |252| , |253|| . Consider a target of external photons 
of energy near E Jjex t with a luminosity L 7iex t- The fraction of proton energy transfered to 
pions is approximately given by 

Aext * UM L T ^l %eXt - (105) 
The neutrino energy threshold is 

VlOOkeVA 



E„ > 7 x 10 13 eV 



J 7,ext 

We will no longer relate the flux of protons to cosmic rays. Instead we introduce the 
luminosity of protons, L p above pion production threshold. This is a largely unknown 
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parameter although it has been estimated to be on the order of 10% of the total luminosity 
|249|| . The proton energy needed to exceed the threshold of Eq. 44 is 

^>1.4xl0-m7«. (107) 



\ 30 / V i? 7l ext 

The neutrino flux can be calculated function of L 



where At is the duration of a blazar flare. This reduces to 

*,~4x 10* km- f { J> , ) f 100 ° MpC V 2 f BY (J^L.) (109) 

V .5 J V10 45 erg/s/V 1000 sec J \ D ) V7/ VlOOkeV/ V ' 

for a fifteen minute flare. Using Eqs. 11, the event rate of TeV neutrinos is 

~ ^ ~ (^) (^)^) 2/5 (^)rno) 

Note that this result is for a typical, but fairly distant source. A nearby line-emitting blazar 
could be a strong candidate for neutrino observation. Considering the more than 60 blazars 
which have been observed, the total flux may generate conservatively tens or, optimistically, 
hundreds of TeV-PeV neutrino events per year in a kilometer scale neutrino telescope such 
as IceCube. The upper range of this estimate can be explored by the AMANDA experiment. 

Blazar neutrino searches should be able to find incontrovertible evidence for cosmic ray 
acceleration in active galaxies, or, alternatively, challenge the possibility that AGN are the 
sources of the highest energy cosmic rays. 



D. Neutrinos Associated With Cosmic Rays of Top-Down Origin 



In addition to astrophysical objects such as GRB and blazars, a variety of top-down 
models have been proposed as the source of the highest energy cosmic rays. For example, 
annihilating or decaying superheavy relic particles could produce the highest energy cosmic 
rays g2§, HH |55|, |56|, |57|, H! |59|, B |6| g62|. In addition to the cosmic ray nucleons, 
they will also generate gamma-rays and neutrinos. Other top-down scenarios which solve the 
cosmic ray problem in a similar way include topological defects ||131| , |132| , |263j , |264| , |265| , |266| 
and Z-bursts |I14|, 115 , [116| , 117 , 118|| . Conventional particle physics implies that ultra 



high-energy jets fragment predominantly into photons with a small admixture of protons 
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267| |268|, [269| . This seems to be in disagreement with mounting evidence that the highest 



energy cosmic rays are not photons |1J. In light of this information, we must assume that 
protons, and not gamma-rays, dominate the highest energy cosmic ray spectrum. This does 
not necessarily rule out superheavy particles as the source of the highest energy cosmic 
rays. The uncertainties associated with the universal radio background and the strength 
of intergalactic magnetic fields leave open the possibility that ultra high-energy photons 
may be depleted from the cosmic ray spectrum near 10 20 eV, leaving a dominant proton 



component at GZK energies [|270| , |271| , |272| , p73j| . With this in mind, one must normalize 
the proton spectrum from top-down scenarios with the observed ultra high-energy cosmic 
ray flux. 

An important point is that this "renormalization" is not only challenged by the large 
sub-GeV photon flux, but also by the neutrino flux associated with these models. Neutri- 
nos, which are produced more numerously than protons and travel much greater distances, 
typically provide an observable signal in operating high-energy neutrino telescopes. 



1. Nucleons in Top- Down Scenarios 

The assumption that nucleons from ultra high-energy fragmentation are the source of the 
highest energy cosmic rays normalizes the rate of their generation. To do this, it is necessary 
to calculate the spectrum of nucleons produced in such jets. Each jet will fragment into a 
large number of hadrons. The quark fragmentation function can be parametrized as |[274|| : 



d ^ = Cx -- i ' 2 {l-xf. (Ill) 
ax 

Here, x = i?Hadron/-Ejet, is the number of hadrons and C = 15/16 is a normalization con- 
stant determined by energy conservation. For a more rigorous treatment of fragmentation, 

sec 



127, 128, 12 



At the energies considered, all flavors of quarks are produced equally. Top quarks im- 
mediately decay into bW ± pairs. Bottom and charm quarks lose energy from hadronization 
before decaying into charmed and other hadrons. Hadrons eventually decay into pions and 
nucleons. The injection spectrum of nucleons produced in parton jets can be approximately 



described by ||275 
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Here, is the number of jets produced per second per cubic meter. N q is the number of 
quarks produced per jet in the energy range concerned, /at ~ .03 is the nucleon fraction in 
the jet from a single quark [|275|1 . 

To solve the ultra high-energy cosmic ray problem, this proton flux must accommo- 
date the events above the GZK cutoff. Observations indicate on the order of 10 -27 events 
m~ 2 s~ 1 sr~ 1 GeV~ 1 in the energy range above the GZK cutoff (5 x 10 19 eV to 2 x 10 20 eV) 
|H, ||. The formalism of a generic top-down scenario is sufficiently flexible to explain the 
data from either the HIRES or AGASA experiments. The distribution of ultra high-energy 
jets can play an important role in the spectra of nucleons near the GZK cutoff. For example, 
the distribution for decaying or annihilating dark matter is likely to be dominated by the 
dark matter within our galaxy. This overdensity strongly degrades the effect of the GZK 
cutoff. 



2. Neutrinos in Top-Down Scenarios 

There are several ways neutrinos can be produced in the fragmentation of ultra high- 
energy jets. First bottom and charm quarks decay semileptonically about 10% of the time. 
Secondly, the cascades of hadrons produce mostly pions. About two thirds of these pions 
will be charged and decay into neutrinos ||275|| . Furthermore, top quarks produced in the 
jets decay nearly 100% of the time to bW^. The W bosons then decay semileptonically 
approximately 10% of the time to each neutrino species. 

Generally, the greatest contribution to the neutrino spectrum is from charged pions. The 
injection spectrum of charged pions is given by ||275|| : 



=< ?S— (H3) 
<J JN 

- 3d -M—^-—. (114) 



The resulting injection of neutrinos is given by |p7i 

r E Jet /N q d£J 

^ 2.34 / -et®^ ( 115 ) 

Using f N ~ .03 and ^f- ~ 1.5 x 10~ 37 , this becomes: 

rEj et /N q ffE EE/ 5 

^(E) ~ 3.0 x 10- 36 / ^ 2 (1 - ^) 2 |ff (116) 

J2.3AE E v Ej et E^ b 
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for each species of neutrino. N q is the number of quarks produced in the fragmentation in 
the energy range of interest. 

To obtain the neutrino flux, we multiply the injection spectrum by the average distance 
traveled by a neutrino and by the rate per volume for hadronic jets which we calculated 
earlier. Neutrinos, not being limited by scattering, travel up to the age of the universe at 
the speed of light (~ 3000 Mpc in an Euclidean approximation). A random cosmological 
distribution of ultra high-energy jets provides an average distance between 2000 and 2500 
Mpc. 

The neutrinos generated in these scenarios can be constrained by measurements of the 
high-energy diffuse flux. AMANDA-B10, with an effective area of ~5,000 square meters has 
placed the strongest limits on the flux at this time. In addition to the diffuse flux of high- 
energy neutrinos, the number of extremely high-energy events can be considered. Depending 
on the details of fragmention and jet distribution, tens to thousands of events per year per 
square kilometer effective area can be generated above an energy threshold of 1 PeV where 
there are no significant backgrounds to interfere with the signal. 

As a simple example, take the Z-burst scenario. In this scenario, ultra high-energy 
neutrinos travel cosmological distances and interact with massive (~eV) cosmic background 
neutrinos at the Z-resonance. The Z bosons then decay producing, among other things, the 
super-GZK cosmic rays. In the center-of-mass frame of the neutrino annihilation the Z is 
produced at rest with all the features of its decay experimentally known. Independent of 
any differences in the calculation, normalizing the cosmic ray flux to the protons, rather 
than the photon flux, raises the sensitivity of neutrino experiments as in all other examples. 
This can be demonstrated with a simple calculation. Data determine that Z-decays produce 
8.7 charged pions for every proton and, therefore, 8.7 x 3 = 26.1 neutrinos from ir^ — > 
v^H — > ev^Vz for every proton. AGASA data, with an integrated proton flux of ~ 5 x 
10 24 ev 2 m~ 2 s _1 sr _1 in the range of, say, 3 x 10 19 to2 x 10 20 eV, indicate ~ .5 protons per 
square kilometer, per year over 2n steridians. This normalization, corrected for the fact 
that neutrinos travel cosmological distances rather than a GZK radius for protons, predicts 
.5 x 26.1 x 3 5Q ( ^p^ c ~ 800 neutrinos per square kilometer, per year over 2ir steridians. We 
here assumed an isotropic distribution of cosmological sources. The probability of detecting 
a neutrino at ~ 10 EeV is ~ .05, see Eq. 11. Therefore, we expect 40 events per year in 
IceCube, or a few events per year in AMANDA II. As with other top-down scenarios, present 
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experiments are near excluding or confirming the model. 

Using the high-energy neutrino diffuse flux measurements and searches for super-PeV 
neutrinos, top-down scenarios can be constrained. Further data from AMANDA, or next 
generation neutrino telescope IceCube, will test the viability of top-down scenarios which 
generate the highest energy cosmic rays. 

IV. THE FUTURE FOR HIGH-ENERGY NEUTRINO ASTRONOMY 

At this time, neutrino astronomy is in its infancy. Two telescopes, one in Lake Baikal and 
another embedded in the South Pole glacier, represent proof of concept that natural water 
and ice can be transformed into large volume Cherenkov detectors. With an acceptance of 
order 0.1 km 2 , the operating AMANDA II telescope represents a first-generation instrument 
with the potential to detect neutrinos from sources beyond the earth's atmosphere and the 
sun. It has been operating for 3 years with 302 OM and for almost 3 years with 677 OM. Only 
1997 data have been published. While looking forward to AMANDA data, construction has 
started on ANTARES, NESTOR and IceCube, with first deployments anticipated in 2002 
and 2003. At super-EeV energies these experiments will be joined by HiRes, Auger and 
RICE. A variety of novel ideas exploiting acoustic and radio detection techniques are under 
investigation, including ANITA for which a proposal has been submitted. Finally, initial 
funding of the R&D efforts towards the construction of a kilometer-scale telescope in the 
Mediterranean has been awarded to the NEMO collaboration. With the pioneering papers 
published nearly half a century ago by Greisen, Reines and Markov, the technology is finally 
in place for neutrino astronomy to become a reality. The neutrino, a particle that is almost 
nothing, may tell us a great deal about the universe [ j277| . 
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